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ABSTRACT 

As  more  and  more  composite  materials  are  used  in 

modern  aircraft  cons t rue t i on ,  the  understanding  of  the 
damage  tolerance  of  this  relatively  stiff,  brittle, 

anisotropic  material  becomes  important  to  designers.  These 
composites  may  suffer  surface  damage  due  to  abrasions  and 
burns,  material  damages  such  as  excess  voids  due  to 

careless  manufacturing  techniques,  or  unseen  damage  in  the 
forms  of  delaminations  due  to  a  low  speed  impact.  These 
damages  all  cause  eccentricities  which  result  in  lower 

panel  buckling  values  when  the  panel  is  compressively 

loaded . 

This  thesis  investigated  the  behavior  of  a  cylin¬ 
drical  composite  panel  made  of  AS4/3501-6  graphite/epoxy 

with  ply  orientations  [O/-45/45/90 3  .  Abrasion  and  burn 

s 

surface  damage  was  physically  modeled  in  the  panels-  by 

removing  a  portion  of  the  exterior  plies.  'The  panels  were 
then  tested  by  compressively  loading  them  and  a  comparison 
was  made  to  buck  ling  predictions  obtained, using  a  STAGSC- 1 
shell  program.  -These  tests  indicated  that  panels  which 

have  suffered  minor  surface  damage  do  not  deviate 

significantly  from  buckling  predictions,  obtained  using  a 
STAGSC-1  linear  bifurcation  model. 


Panels  were  also  tested  which  had  varying  thicknesses 


* 


and  variations  in  void  content  due  to  faulty  manufacturing 
techniques.  These  panels  were  also  compres s i ve 1 y  loaded, 
and  it  was  found  that  high  void  content  increased  panel 
thickness,  which  resulted  in  higher  compressive  strengths. 

'Composite  laminates  subjected  to  a  low  speed  impact, 
such  as  a  dropped  tool  or  a  manufacturing  load,  often 
develop  an  internal  de lamina t i on .  This  delamination  may 
result  in  the  reduction  of  the  panel’s  strength  when 
subjected  to  compressive  load. 

'“’Since  curved  panels  are  3-dimensional,  and  buckling  is 
a  non-linear  phenomenon,  the  compressive  load  which  will 
cause  curved  panels  to  become  unstable  is  extremely  hard  to 
predict  analytically.  This  thesis  presents  a  technique 
whereby  the  local  buckling  loads  at  the  delamination  may  be 
predicted  using  a  2-dimens ional  model  with  a  plane  strain 


correction.,.  This  model  yielded  predictions  of  local 


instability  within  30  percent  of  experimental  values. 
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A  STUDY  OF  DAMAGE  TOLERANCE 
IN  CURVED  COMPOSITE  PANELS 

I .  Introduction 

Composite  materials  are  being  used  more  and  more 
frequently  in  the  design  and  construction  of  modern 
aircraft.  Composites  are  used  for  aircraft  skins,  for  major 
structural  components  and  for  small  component  parts. 

Composites  have  many  advantages  over  the  more 

homogeneous,  isotropic  metals  such  as  aluminum  and  titanium 

that  have  been  used  in  classical  aircraft  construction. 

Modern  designers  can  choose  from  a  wide  variety  of 

composite  fibers  and  matrix  resins  to  tailor  their  strength 

and  stiffness  to  a  given  application.  Composites  can  be 

chosen  which  have  higher  strength  at  high  temperatures  than 

metals,  or  which  are  more  dimensionally  stable  over  a  wider 

temperature  range  than  metals.  Composites  have  a  much 

higher  specific  stiffness  and  specific  strength  than 

6  6 

metals,  with  stiffness  ranging  from  2x10  psi  to  18x10 
psi  and  strengths  ranging  from  10  ksi  to  over  200  ksi. 
Fatigue  limits  are  far  in  excess  of  those  in  aluminum,  with 
a  much  greater  vibrational  damping. 

The  ability  to  achieve  great  strength  at  reduced 
weight  has  had  a  significant  impact  on  spacecraft  and 
aircraft  design.  Composites  are  currently  being  used  in 


a i rc  raft 


such  as 


the  Grumman  X-29  Advanced  Technology 


Fighter,  the  General  Dynamics  F-16  Falcon  and  the  McDonnel 
Douglas  F-15  Eagle.  They  are  also  being  used  in  the  Spare 
program  where  weight  savings  and  dimensional  stability  over 
a  wide  temperate^’  range  are  critical.  Currently,  the 
world’s  largest  composite  structures  are  the  carbon-fiber 
reinforced  motor  casings  for  the  space  shuttle’s  solid  fuel 
rockets.  Composites  are  also  used  in  space  frames  such  as 
the  shuttle’s  retractable  boom,  and  in  the  trusses  used  to 
stabilize  the  spot  beam  reflectors  for  the  Intelesat  V 
communications  satellite.  [5] 

Graphite  epoxy  laminates  are  perhaps  the  most 
versatile  and  prevalent  of  composite  materials.  They  have 
been  in  wide  use  for  a  number  of  years  and  are  readily 
available.  In  the  early  1970’s,  the  increased  use  of 
graphite  epoxy  in  commercial  products  such  as  fishing  rods, 
skis,  and  tennis  racquet  has  resulted  in  a  continuing  drop 
in  the  cost  of  graphite/epoxy.  Not  only  are  they  becoming 
less  expensive,  but  they  can  give  strength  comparable  to 
aluminum  at  weight  savings  of  15  to  30  percent.  With  the 
use  of  computer  controlled  and  advanced  robotics,  automated 
filament  windings  are  possible.  As  a  result,  parts  can  be 
consolidated,  and  the  post  fabrication  machining  often 
used  for  metal  parts  can  be  reduced.  Up  to  50%  reductions 
in  fabrication  costs  are  the  end  result. 

However,  there  are  also  potential  disadvantages  to  the 
use  of  composites.  The  design  of  composite  structures  is 


much  more  complicated,  and  less  intuitive  than  design  with 
isotropic  materials.  Careful  consideration  must  be  given 
to  the  rotational  effects  caused  by  variations  in 
temperature  and  ply  lay-up  in  a  composite.  A  composite 
panel  can  peel  due  to  air  friction  in  a  high  performance 
aircraft,  and  damage  to  composite  structures  is  frequently 
more  complicated  and  time  consuming  to  repair  than  damage 
to  metal  structures.  Damage  to  a  composite  caused  by  a  low 
velocity  impact  such  as  a  dropped  tool,  manufacturing  loads 
or  strikes  from  runway  debris  is  not  easily  detected  and 
can  cause  serious  reduc  tions  in  the  strength  of  the 
structure  resulting  in  the  collapse  of  the  structure  at 
loads  which  are  much  smaller  than  the  design  load. 

BACKGROUND 

The  problem  of  shell  collapse  has  been  discussed  in 
the  literature  for  many  years,  and  analytic  techniques  have 
been  developed  to  predict  the  buckling  load  for  simple 
isotropic  materials  with  simple  geometries.  Stability 
equations  for  cylindrical  shells  were  first  developed  in 
the  late  1800’s.  In  1911,  a  closed  form  solution  for  the 
collapse  of  a  cylindrical  shell  was  presented  by  Lorenz. 
In  1932,  Flugge  presented  a  comprehensive  treatment  of 
cylindrical  shell  stability  for  shells  subjected  to 
compression  and  bending  loads  [2].  Fl.gge  tested  a  series 


of  shells  and  found  that  his  experimental  values  were 
approximately  50  percent  of  his  closed  form  predictions. 
In  1933,  Donnell  presented  a  study  which  treated  the  class 
of  shells  known  as  quasi-shallow  shells.  These  shells  were 
relatively  flat  as  are  most  structural  shells,  and  because 
of  the  simplicity  of  his  solution,  Donnell’s  equat i ons 
became  widely  used.  Donnell’s  experimental  values  for  thin 
cylindrical  aluminum  shells,  however,  were  even  farther 
from  his  predictions  than  Fl.gge's  had  been  [2].  Finally, 
von  Karman  and  Tsien  performed  an  analysis  in  1941  on  the 
post  buckling  paths  of  compressed  cylindrical  shells. 
Their  work  showed  that  the  secondary  path  dropped  sharply 
downward  from  the  bifurcation  point.  As  a  result,  slight 
imperfections  in  a  shell  induced  rotational  effects  which 
caused  the  actual  equilibrium  path  to  fall  at  a  lower  load 
level  than  the  theoretical  linear  path.  Hence,  for  shells 
with  a  slight  imperfection,  the  actual  collapse  values 
could  be  quite  small  compared  to  predictions. 

As  composite  panels  came  on  the  scene,  the  problem  of 
shell  instability  was  compounded.  A  composite  panel  which 
contained  a  delamination  often  had  a  thin  section  of 
delaminated  plies  attached  to  a  thicker  base  laminate.  As 
a  result,  the  delaminated  plies  were  more  flexible  than  the 
base  laminate,  and  were  often  the  first  portion  of  the 
laminate  to  experience  instability. 


A 


Early  work  with  de larainat ions  in  composite  materials 
was  done  by  attaching  a  thin  composite  layer  to  a  thick 
isotropic  base  material  such  as  aluminum.  This  technique 
allowed  the  introduction  of  pure  compression  loads  without 
having  to  worry  about  rotational  effects  in  the  base 
material.  In  1981,  Whitcomb  [17]  used  this  technique  to 
study  a  rectangular  delamination  in  a  graphite/epoxy 
panel.  He  used  a  thin  teflon  film  to  prevent  bonding  of 
the  delaminated  region  to  the  base  material.  Using  this 
technique,  he  was  able  to  study  Mode  I  and  Mode  II  strain 
energy  release  rates  and  their  effects  on  interlaminar 
stresses  at  the  delamination  crack  tip. 

In  1984,  Whitcomb  [18]  substituted  the  aluminum  base 
material  for  a  thick  quasi - isotropic  graphite/epoxy  base 
laminate.  Still,  using  a  piece  of  teflon  to  form 
rectangular  strip  delaminations  in  a  flat  plate,  Whitcomb 
found  that  buckling  can  cause  high  interlaminar  stresses  at 
the  delamination  crack  tip  and  that  buckling  strain  is  an 
important  phenomenon  in  assessing  potential  for  strength 
loss  due  to  a  delamination.  In  this  study,  Whitcomb  [18] 
analyzed  the  complicated  non-linear  buckling  phenomenon 
using  the  superposition  of  several  minor  linear  problems. 

Actual  delaminations  caused  by  an  impact  were  found  to 
be  elliptical  in  shape,  however.  These  elliptical 
delanunations  in  flat  plates  were  studied  by  Shivakumar  and 
Whitcomb  in  1984  [13].  This  study  also  used  flat  quasi¬ 
isotropic  plates,  and  developed  a  simple  Rayleigh  Ritz 


solution  based  on  the  Trefftz  criterion  to  predict  local 
buckling  of  delaminated  plies  (referred  to  as  the 
sublaminate)  from  a  quas i - i so t ropic  base  laminate.  They 
found  that  there  were  no  interlaminar  stresses  in  the 
region  of  the  delamination  until  the  onset  of  buckling  and 
that  the  size  of  the  delamination  was  an  important 
parameter  in  the  stress  intensity  at  the  de 1  am  1 na t i on . 
However,  this  technique  gave  inaccurate  results  for  highly 
anisotropic  sublaminates,  and  the  technique  wasn’t  extended 
to  curved  shells.  This  study  attempts  to  address  the 
problem  of  the  buckling  of  an  anisotropic  ( 0 /  —  4  5 } 
sublaminate  in  a  curved  composite  shell. 

All  of  the  studies  involving  composite  materials  that, 
were  mentioned  above  had  involved  a  flat  plate  geometry. 
Since  most  real  structural  shells  are  not  flat, 
generalizations  made  from  these  studies  sometimes  couldn’t 
be  applied  to  a  curved  surface.  In  the  early  1970's,  a 
researcher  at  General  Dynamics  named  Wilkens  began  using  a 
fixture  to  test  curved  laminate  panels  [20],  About  this 
time,  Dr.  Anthony  Palazotto  initiated  research  at  AFIT  into 
the  behavior  of  curved  laminate  panels.  In  1977,  Nelson 
studied  the  buckling  behavior  of  curved,  stiffened  panels 
with  cutouts.  In  1979,  Becker  compared  experimental 
buckling  loads  in  composite  panels  to  analytic 
predictions.  Later,  Bauld  investigated  the  effects  of  the 
aspect  ratio  of  the  panels  on  buckling  loads.  In  1981, 


another  AFIT  student,  Siefert  investigated  the  effects  of 
centrally  located  delamimtions  on  composite  panel  strength 
[12].  Horban  [7]  continued  this  research  in  1985  by  using 
a  non-linear  STAGS  finite  element  code  to  study  multiple 
de lam l nat i ons  in  curved  panels. 

Composite  panels  present  other  unique  problems  to  the 
designer.  Laminates  tend  to  develop  large  interlaminar 
stresses  which  depend  to  a  large  extent  on  the  orientations 
of  the  laminate’s  fibers  and  the  locations  of  singularities 
such  as  holes  or  boundaries.  Due  to  the  expense  and 
complexity  of  computer  models  required  to  model  a 
delamination  in  a  composite  panel,  most  of  the  work 
invoicing  discontinuities  in  a  composite  panel  have  been 
done  analytically.  However,  analytic  results  have  only 
been  developed  for  a  narrow  range  of  geometries  which 
neglect  any  defects  in  the  material  such  as  surface 
imperfections  or  internal  de  1  am  l  na  t  i  ons  .  This  study- 
attempts  to  address  this  problem  by  developing  a  technique 
whereby  the  effects  of  a  delamination  in  a  curved  composite 
panel  can  be  modeled  using  relatively  inexpensive  linear 
shell  models  with  a  plane  strain  correction. 


PURPOSE 


The  purpose  of  this  thesis  was  to  experimentally  study 
the  tolerance  of  a  curved  composite  panel  to  damage. 
Damages  investigated  included  surfac  damage  due  to 
abrasion  or  burning,  higher  than  normal  void  contents  due 
to  faulty  manufacturing,  and  delamination  damage  due  to  low- 
speed  impact.  This  study  also  attempted  to  develop  an 
inexpensive  two  dimensional  model  using  the  superposition 
of  linear  buckling  finite  element  codes  to  predict  the  non¬ 
linear-  three  dimensional  phenomenon  of  sublaminate  buckling 
in  a  thin  curved  panel.  This  investigation  involved  both 
analytical  work  and  experimental  work  to  verify  the 
predicted  local  buckling  values. 

Internal  de lamina t i ons  in  a  composite  panel  can 
develop  as  a  result  of  a  low  velocity  impact  caused  by  a 
dropped  tool,  a  projectile,  or  as  a  result  of  manufacturing 
loads.  These  delaminations  weaken  the  panel  and  reduce  the 
load  at  which  panel  "global”  buckling  will  occur.  Another 
phenomenon  that  has  been  noted  [7]  is  the  snapping  effect 
of  the  delaminated  plies  as  they  reach  their  local  buckling 
load.  This  snap  out  may  occur  at  a  load  which  is  much  less 
than  the  global  buckling  load  of  the  panel.  Since  curved 
composite*  panels  are  3 -d  i  mens  i  ona  1  and  buckling  is  a  non- 
l.near  problem,  the  compressive  load  which  will  cause  the 
panel  to  become  unstable  is  extremely  hard  to  predict, 
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analytically.  Another  purpose  of  this  study  is  to  develop 
a  method  which  will  use  a  2 -d i mens i ona 1  linear  finite 
element  model  and  a  plane  strain  correction  to  accurately 
predict  this  initial  instability. 

SCOPE 

The  panels  used  in  this  study  were  quasi-isotropic 
cylindrical  composite  panels  made  of  AS4/3501-6 
graphite/epoxy  with  ply  orientations  of  ( 0/ -4 5 /  +  4 5 /90  )  . 
Their  height  was  12  inches  and  they  had  a  radius  of 
curvature  of  12  inches,  (fig  1.1)  The  actual  experimental 
panels,  however,  were  trimmed  slightly  oversize  to  allow 
for  mounting  in  a  test  fixture. 

Previous  work  done  by  Horban  and  Palazotto  [7]  has 
shown  that  a  teflon  film  placed  between  plies  during  panel 
layup  will  cause  a  total  internal  delamination  in  the 
panel.  Therefore,  0.5  mil  teflon  discs  were  placed  at  the 
geometric  centers  of  the  panels  between  plies  2-3  on  the 
panels’  concave  side,  and  between  plies  6-7  on  the  convex 
side  (fig  1.2).  Panels  were  manufactured  with  a  four  inch 
and  two  inch  diameter  delamination  and  with  no  dc lami¬ 
na  t  i  oris  . 

Panels  were  also  manufactured  with  a  teflon  disc 
placed  between  plies  1-2  and  between  plies  7-8.  The 
delam i na t ions  caused  by  placing  inserts  adjacent  to  the 


exterior  pLy  resulted  in  delaminations  that  were  externally 
visible.  The  delamination  appeared  to  be  a  circular 
blister  approximately  the  size  of  the  insert,  with 
fractures  in  the  matrix  running  parallel  to  the  fiber 
orientation  of  the  blistered  ply.  Since  the  blistered 
portion  of  the  external  ply  could  take  no  compressive  load, 
it  was  removed  and  the  panel  was  used  to  model  a  panel 
which  has  suffered  external  damage  resulting  from  abrasion, 
chafing,  or  burning.  These  panels  were  then  experimentally 
tested  to  determine  the  effect  of  this  damage  on  panel 
strength . 

The  boundary  conditions  (see  fig  1.1)  used  in  the 
experimental  fixture  were:  simply  supported  sides 


(u=v=w  =free  w=w^=0),  a  clamped  bottom  edge  ( u=v=w=w^=0 ) , 
and  a  clamped  top  edge  (u=free  v=w=w  =0).  A  constant 
displacement  load  was  introduced  at  the  top  edge  of  the 


panel  . 


1 


Concave 
Inside  Surface 


Figure  1.2  -  Cross  Section  of  Panel 
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II.  THEORY 


Classical  Laminated  Plate  Theory 


An  overview  of  classical  laminate  plate  theory  is 
presented  here  for  readers  not  familiar  with  the  behavior  of 
composite  materials  and  the  associated  notation.  Interested 
readers  may  read  Jone’s  text  [8]  for  a  more  in-depth 
treatment  of  this  material.  This  section  is  a  condensation 
of  Jone’s  text,  Section  II  of  the  SQ5  user  manual  [11]  and 
some  work  presented  in  previous  theses  [6,7,12]. 

Co-ordinates  in  an  orthotropic  layer  of  a  laminate  are 
defined  with  respect  to  the  material's  fiber  orientation. 
The  1-direction  is  defined  as  the  direction  parallel  to  the 
fiber,  and  the  2-direction  is  the  direction  normal  to  the 
fibers  in  the  plane  of  the  laminate  layer  (see  fig  2.1). 

To  form  the  constitutive  relations  in  terms  of  the  1-2 
directions,  the  four  independent  material  properties 
( ,E2 ,G12 , )  are  used: 


Q1 1  Q12 
Q12  Q22 
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Figure  2.2  -  Geometry  of  N  Layered  Laminate 
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When  the  designer  is  interested  in  stresses  in 

directions  other  than  the  1-2  directions,  then  the  Q.  .  matrix 

i  J 

is  transformed  by  the  angle  <9  between  the  1  axis  and  the  new 

x  axis  of  interest  (see  fig  2.1).  When  this  transformation 

is  performed  for  each  Q.  .  term,  the  following  relation 

^  J 

results: 
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(2.3) 


Where  the  Q  are  the  stiffness  terms  transformed  by  the 
rotation  0  of  the  1-2  axes  into  the  x-y  axes. 

The  displacements  at  any  point  in  the  laminate  cross 
section  are  defined  in  terms  of  the  strains  at  the  midsurface 


o  o  o 


of  the  laminate  (Key)  and  the  plate  curvatures  (  k  k 

x  y  xy  x  y 
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(2.4) 


For  the  relation  given  in  (2.4),  a  coupile  of  assumptions 
must  be  made.  First,  the  laminae  are  assumed  t.o  be  very  thin 
with  displacements  continuous  across  the  laminae  thickness. 
Second,  the  Ki rchho f f-Love  assumption  applies.  This 
assumption  says  that  normals  to  the  mid-surface  remain  plane 
and  normal  to  the  mid-surface  after  bending.  For  relatively 
thin  laminates  with  small  rotations  relative  to  the  neutral 
axis,  this  is  not  a  bad  assumption.  For  a  cylindrical  shell, 
the  strains  and  curvatures  expressed  in  equation  (2.4)  above 
are  defined  as  follows,  with  strains  at  the  midsurface  given 
by  : 
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and  curvatures  given  by: 


k  =  0  x  ,  x 
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kxy  =  2  L  0y’X  +  0X,Y  +  R 


(2.6) 


In  equations  (2.5)  and  (2.6)  above,  R  is  the  panel’s 
radius  of  curvature;  u,v,w  are  the  axial,  circumferential, 


and  radial  displacements  of  the  shells  mid-surface;  and  the 


t 


L 


3  9  s  are  rotational  components  which  are  expressed  in  terms  of 


displacement  as  follows: 
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(2.7) 


1  /  2  v  >  x  -  u  ,  y 


Now  that  the  stresses  have  been  computed  for  each  individual 

ply  using  equation  (2.3),  the  force  resultants  (N  ,N  ,N  ), 

x  y  x  y 

and  the  moment  resultants  (M  MM)  can  be  computed  by 

x  y  xy 

integrating  the  laminae  stresses  over  the  laminate  thickness 
in  the  z  direction. 


(2.8) 


z  dz 


(2.9) 


Now  we  substitute  equation  (2.3)  into  (2.8)  and  (2.9). 
In  matrix  notation,  we  get  the  following  equations: 
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Defining  the  quantities: 
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We  can  now  write  equations  (2.10)  and  (2.11)  in  terms  of 
the  extensionai  stiffness  matrix  [A],  the  coupling  stiffness 
matrix  [B],  and  the  bending  stiffness  matrix  ID),  where  IA], 
(B)  and  [D]  are  each  a  3x3  matrix  whose  values  are  given  in 
appendix  A.  Now  we  can  write  an  expression  for  N  and  M  in 
matrix  notation: 


(2.15) 


Equatioin  (2.15)  is  the  total  constitutive  relation  for 
a  laminated  plate.  The  coupling  of  inplane  extension  and 
bending  is  seen  to  be  a  result  of  the  [B)  matrix.  lor  a 
laminate  that  is  symmetric  about  the  midplane,  the  |B  ] 
terms  are  zero,  and  there  is  no  coupling  of  extension  with 
bend i ng . 
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a  structure  wi  11  reach  when  second  order  rot  ational  effect  s 
are  neglected  and  is  the  point  at  which  the  structure  will 
become  unstable  and  buckle. 

Figure  2.4  represents  the  equilibrium  path  for  a 
perfect  circular  cylindrical  shell.  No  structure  is 
perfect,  however,  and  for  the  case  of  cylindrical  shells, 
very  slight  shell  imperfections  result  in  large  reductions 
of  the  buckling  load.  Experimental  collapse  loads  for 
axially  compressed  cylindrical  shells  may  even  be  as  low  as 
20%  of  the  theoretical  values  [1]. 

The  reason  for  this  large  difference  in  theoretical  vs 
actual  load  can  be  seen  by  observing  figure  2.4.  For 
circular  cylindrical  shells,  the  secondary  equilibrium  path 
drops  sharply  downward  from  the  bifurcation  point.  Since 
the  equilibrium  path  is  the  path  for  a  perfect  structure,  a 
slight  imperfection  in  the  shell  will  induce  eccentricities 
which  increase  the  second  order  rotational  effects  and 
result  in  the  actual  curve  falling  somewhat  below  the  ideal 
curve.  Therefore,  the  critical  load  will  occur  at  a 
smaller  value,  and  no  post  buckling  increase  in  load  can 
result.  The  reduction  in  the  actual  strength  of  a  shell 
from  the  critical  load  at  bifurcation  is  known  as  knockdown 
and  is  essentially  the  percent  difference  between  predicted 
values  and  experimental  values.  Horban  [7]  and  Siefert 
[12]  used  this  knockdown  factor  in  their  work  with 
cylindrical  composite  panels  to  compare  experimental 
results  with  the  higher  STAGSC-1  bifurcation  predictions. 


Load 
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Figure  2.4  -  Equilibrium  Path  for  Cy I i nderical  Shell 
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STAGSC-1  Theory 


STAGSC-l  (Structural  Analysis  of  General  Shells)  is  an 
energy  based  finite  element  program  for  thin  shells. 
STAGSC-l  began  its  development  in  1967  under  government 
contract  to  Lockheed  La  1 o  Alto  Research  Laboratory,  and  its 
development  has  continued  for  the  past  twenty  years  under 
various  contracts  through  NASA  Langley  Research  Center. 
The  first  version,  known  as  STAGS,  was  a  finite  difference 
program  for  the  non-linear  analysis  of  thin  shells  with 
u touts.  A  linear  version  for  shells  of  revolution  was 

developed  in  1968,  and  buckling  and  thermal  capabilities 
were  added  in  1970.  The  code  evolved  further  in  1970  with 
the  add i t ion  of  more  elements  to  its  library  and  an 
inelastic  capability.  In  1973,  a  version  called  STAGSA  had 
evolved  which  included  transient  response,  dynamic 

eigenvalue  analysis,  and  dynamic  buckling  response 

routines.  STAGSC-l  was  released  in  1979.  The  new  code  was 
entirely  transformed  to  a  finite  element  based  program,  and 
was  particularly  useful  for  bifurcation  analysis  and 
dynamic  analysis  of  thin  shells  [  1  4  j  . 

The  lutes'  version  of  STAGSC-l  was  released  in  1986 
and  addresses  the  full  collapse  of  a  structure.  Also, 
previous  versions  of  STAGS  were  limited  to  problems  in  the 
small  rot  at  ion  class  (less  than  10  degrees).  Ihe  present 
version  addresses  this  problem  using  an  element -i ndependent 


co-rotational  formulation  for  ail  static  non -linear 
collapse  problems,  and  it  can  now  be  used  for  structures 
with  large  rotations  (greater  than  10  degrees)  at  small 
strains.  The  1986  version  also  uses  the  Riks  algorithm  for 
non- linear  collapse  evaluation  typical  of  optimized  shells 
made  of  the  new  lightweight  materials  (10).  The  new  Riks 
method  is  an  efficient  nonlinear  solver  which  can  extend 
STAGS  solutions  into  the  post  buckled  portion  of  the  load 
versus  deflection  curve.  The  iy86  program  is  the  version 
that  is  used  for  this  research. 

The  ST.AGSC-1  computer  program  is  made  up  of  four  main 
parts:  STAGS1,  STAGS  2 ,  POSTP,  and  STAPL.  The  STAGS1 
portion  is  the  part  of  the  program  which  compiles  the  input 
deck,  performs  model  generation,  and  begins  performing 
preliminary  calculations.  The  STAGS2  module  does  the  bulk 
of  the  numerical  calculations  performed  in  a  finite  element 
analysis.  STAGS2  performs  such  tasks  as  matrix 
decomposition,  linear  and  non-linear  stress  analysis, 
eigenvalue  analysis,  thermal  gradient  determination,  and 
the  calculation  of  stress  resultants,  moments  and 
deflections.  POSTP  is  a  post-processor  for  determining 
secondary  solutions  from  previously  calculated 
displacements  which  have  been  saved  to  TAPE22.  The  STAPL 
module  is  a  plotting  routine  which  can  be  used  to  plot  the 
undeformed  finite  element  mesh,  the  deformed  mesh,  and 
contour  plots  of  the  structure’s  deflections  (19(.  Readers 
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interested  in  running  the  STAPL  module  should  refer  to  the 
STAGSC-1  user  manual  [1|  and  to  Tisler  [15),  appendix  B  for 
sample  input  decks  for  restart  and  STAPL  runs. 

The  1986  version  of  STAGSC-1  that  was  used  for  this 
research  has  two  routines  which  are  both  called  STAGS2. 
Either  routine  can  be  loaded  in  the  computer  under  the  name 
STAGS2.  The  first  routine  performs  a  linear  bifurcation 
analysis  state,  while  the  second  performs  its  computations 
using  a  non-linear  buckling  routine.  "It  is  suggested  that 
users  of  STAGSC-1  utilize  the  non-linear  module,  as  it  will 
give  more  accurate  results  at  approximately  the  same  cost 
in  computer  time.  True  linear  problems  will  converge  very 
quickly  using  the  non-linear  version,  so  there  is  usually 
no  real  advantage  in  running  a  linear  model."  [10J  If  a 
user  should  encounter  problems  with  a  model  that  does  not 
run,  he  should  check  the  B1  card  in  the  input  deck  to 
ensure  that  the  analysis  being  performed  i s  compatible  with 
the  STAGS2  module  called  on  the  system  [19]. 

Another  unique  aspect  of  STAGSC-1  is  its  use  of 
surface  co-ordinates  on  the  shell  surface.  This  allows  any 
point  on  the  mesh  to  be  defined  in  terms  of  two  independent 
surface  co-ordinates;  or  conversely,  a  user  can  use  the 
system’s  global  frame  to  define  the  grid  points.  Also, 
shell  units  consisting  of  various  elements  and  material 
properties  may  be  formed.  These  shell  units  can  then  be 
interconnected  using  the  STAGS  G1  card  to  form  a  continuous 
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compatible  shelL.  This  feature  is  useful  when  modeling 
stiffeners,  beams,  or  shells  whose  geometries  change  along 
the  surface . 

The  STAGSC - 1  program  calculates  the  bifurcation  point 
using  variational  potential  energy  techmaMes .  The  static 
equilibrium  path  is  calculated  using  minimum  potential 
energy  such  that  static  equilibrium  is  achieved  for  a  given 
input  load  or  deflection  (summation  of  external  forces 
equal  zero  at  equilibrium).  Potential  energy,  V,  is  the 
internal  strain  energy,  1,  minus  the  work  done  by  external 
forces,  W;  where  work  is  the  product  of  external  forces  and 
their  associated  deformations. 

V  =  U  -  W  (2.16) 

U  =  1  /  2  {  )  T  [  N  ]  {  -  )  dA  (2.17) 

w  =  ^ { x ) T { F )  (2.18) 

where:  {  k  )  =  vector  of  strains  and  curvatures 

[NJ  =  matrix  of  A.  B and  D.  .  terms 

x  J  i  J  i  J 

{ x )  ^  vector  of  displacements 

(F)  =  vector  of  applied  external  forces 

For  static  equilibrium,  the  potential  energy,  V,  must 
be  a  minimum;  therefore,  the  first  variation  of  potential 


energy,  ;V,  must  equal  zero.  At  the  bifurcation  points  foi 
a  conservative  mechanical  system,  the  second  variation  of 
potential  energy,  f*V  must  be  equal  to  zero.  The  :-V  terms 
represent  the  second  order  potential  energy  terms  which 
result  from  a  small  virtual  displacement  to  the  structure. 
The  critical  load  for  a  structure  occurs  when  the  structure 
first  loses  its  stability.  For  stability  in  a  structure, 
the  second  variation  must  be  positive  definite. 

STAGSC-1  performs  an  eigenvalue  solution  to  find  the 
bifurcation  point.  Readers  are  referred  to  Brush  and 
Almroth  [2]  for  an  in-depth  development  of  potential  energy 
formulations  and  bifurcation  analysis,  and  to  Rankin  [10] 
for  the  1986  updates  to  the  STAGSC-1  computer  code. 


III.  MANUFACTURING  k  EXPERIMENTAL  TECHNIQUES 


Introduction 


Forty  curved  panels  were  fabricated  and  tested  for  this 


research.  The  panels  were  constructed  of  Hercules  AS4/3501-6 


Graph i te / Epoxy  12  inch  prepreg  tape.  All  panels  were  laid 


up  with  a  ( 0  / - 4 5 / 4 5 / 90 )  quasi-isotropic  geometry.  The 


panels  were  cylindrical  shells  having  a  radius  of  curvature 


of  12  inches  measured  to  the  outside  convex  surface  of  the 


panel.  They  had  a  trimmed  height  of  13  inches,  and  an  arc 


length  of  12.5  inches  (fig  1.1).  The  material  properties 


were  assumed  to  be: 


=  18.844  x  10  psi 


1.468  x  10  psi 


=  0.91  x  10  psi 


0 . 28 


=  0.022 


Four  of  the  panels  were  built  with  no  internal 


de 1 am i nat i ons .  These  panels  were  to  serve  as  control  panels 


to  check  the  accuracy  of  computer  models,  and  to  compare 


with  previous  work  done  by  Horban  f  7 )  and  Siefert  [12] 


Internal  de 1  am i na t i ons  caused  by  low  speed  impacts  were 


modeled  by  incorporating  0.5  mil  teflon  discs  between  plies 


during  lay  up.  Two  inch  diameter  and  four  inch  diameter 
delaminations  were  used. 

Previous  work  [7, 12]  had  shown  that  a  complete 
delamination  will  occur  in  a  panel  if  a  teflon  insert  is 
placed  in  the  panel  during  lay  up.  This  work  indicated  that 
it  was  desirable  to  use  a  thin  layer  of  material  (0.5  mil 
mylar)  to  avoid  matrix  fracture  when  the  deiaminations  are 
placed  within  one  ply  of  the  outer  surface  of  the  panel. 
Horban  also  found  that  mylar  inserts  only  resulted  in  a 
partial  de  1  am  i  na  1. 1  on  ,  while  teflon  discs  resulted  in  a  total 
delamination  in  the  panel.  Therefore,  0.5  mil  teflon  was 
chosen  as  the  delamination  causing  insert  for  this 
research . 

Initially,  delaminations  were  caused  by  placing  two 
discs  of  DuPont  300  A,  0.5  mil  FEP  teflon  film  back  to  back 
between  plies  during  lay  up.  The  discs  were  coated  with  RAM 
225  release  agent  prior  to  being  placed  in  the  panel  to 
ensure  100%  debonding  .  Originally,  delaminations  were  to 
be  placed  between  plies  1-2,  2-3,  6-7  and  7-8  (fig  1.2). 
The  inserts  placed  between  plies  2-3,  and  6-7  formed  the 
desired  delamination.  However,  problems  were  encountered 
when  the  teflon  discs  were  placed  only  one  ply  in  from  the 
exterior  surface  (between  plies  1-2  and  7-8).  After  the 
panels  were  removed  from  the  autoclave,  the  deiaminations 
were  externally  visible  (fig  3.1).  They  appeared  to  be  a 
low,  circular  blister  approximately  the  size  of  the  teflon 


insert,  with  fractures  in  the  matrix  running  parallel  to  the 
0  fiber  orientation  of  the  exterior  blistered  ply. 

Horban  [7J  had  previously  noted  this  problem  when  he 
used  one  mil  mylar  inserts.  He  had  then  tried  using  0.5  mil 
mylar  inserts,  and  found  that  he  observed  no  matrix  cracking 
and  no  blistering  effect.  As  a  result  of  Horban’s 
experience,  it  was  thought  that  using  one  0.5  mil  teflon 
insert  with  RAM  225  release  agent  as  opposed  to  the  two 
inserts  placed  back  to  back  would  give  acceptable  results. 

Four  panels  were  manufactured  in  this  manner.  Two  of 
the  panels  had  two  inch  diameter  inserts  while  the  other  two 
had  four  inch  diameter  inserts.  All  four  panels  exhibited 
the  matrix  cracking  and  blistering  seen  previously.  At  this 
point  it  was  not  known  what  was  causing  the  blisters  to 
form.  Panels  were  constructed  using  one  0.5  mil  teflon  disc 
with  no  RAM  225  release  agent  applied.  Blisters  still 
occurred.  It  was  then  thought  that  perhaps  volatiles  in  a 
gaseous  state  were  being  trapped  beneath  the  teflon  disc 
resulting  in  a  vapor  pressure  as  the  panels  were  heated  in 
the  autoclave.  Panels  were  then  constructed  using  one  0.5 
rail  teflon  disc  that  had  been  perforated  with  20  small 
holes.  Blisters  still  occurred.  A  finite  element  model  was 
run  to  determine  if  residual  curing  sresses  might  be  causing 
this  blistering.  No  residual  stresses  were  found  in  the 
model.  Due  to  time  constraints,  and  the  number  of  panels 
a  1  lowed  for  this  research,  no  further  experimental 
investigation  was  attempted.  However,  it.  was  decided  to  u:-  e 


the  panels  with  the  blistered  de  1  am i na t i on  to  carry  out 
further  experimental  work.  Since  the  blistered  portion  of 
the  external  ply  could  take  no  load,  the  blistered  portion 
was  removed  using  a  razor  knife,  and  the  panel  was  used  to 
model  a  panel  which  has  suffered  external  damage  as  a  result 
of  abrasion,  chafing  or  burning.  These  panels  were  then 
experimentally  tested  to  determine  the  effect  of  this  damage 
on  the  panel's  buckling  strength. 

Another  type  of  damage  studied  in  this  thesis  is  panel 
damage  due  to  manufacturing  var i ab i 1 i t i es .  Panels  were 
tested  that  had  laminate  thicknesses  ranging  from  39.3  to 
42.6,  and  which  had  void  contents  ranging  from  0.10  percent 
to  4.4  percent.  These  variations  could  result  from  faulty 
manufacturing  procedures  or  from  mechanical  difficulties 
such  as  a  faulty  vacuum  during  curing  in  the  autoclave. 

Panel  Lay  Up  and  Curing 

The  panels  were  laid  up  in  curved  steel  molds  with  a 
12  inch  radius  of  curvature  (fig  3.2)  .  The  12  inch  prepreg 
was  cut  and  placed  as  shown  in  the  fabrication  plan  (figs 
3. 5-3.7).  Panels  that  were  to  receive  teflon  inserts  were 
carefully  measured  to  the  center  of  the  intended 
de lami nation  using  a  steel  scale,  and  the  teflon  was  placed 
using  a  RAM  225  release  agent  to  ensure  total  de 1  am i na 1 1 on . 


AFIT  Curved  Panels 

15.5"  X  32"  Panel  -  12"  wide  mat’! 

0°  Plies: 


For  One  Panel: 

Cut  2  -  40"  long  X  12"  wide  pieces 
Cut  2  -  40"  long  X  8"  wide  pieces 
(Save  the  extra  4"  wide  X  40"  long  piece  for  the  90°  plies) 

For  All  17  Panels: 

Cut  34  -  40"  long  X  12"  wide  pieces 

Cut  34  -  40"  long  X  8"  wide  pieces  (Save  the  4"  wide  X  40"  lo 
pieces  for  the  90°  plies) 


Figure  3.4  -  Prepreg  Cutting  Plan  for  0  Plies 


AFIT  Curved  Panels 

15.5"  X  32"  Panel  -  12"  wide  mat'l 

90°  Plies: 


For  One  Panel: 

Cut  6  -  20"  long  X  12"  wide  pieces 

Trim  the  4"  wide  X  40"  long  piece  into  2  -  4"  wide  X  20"  long  pieces 
For  All  17  Panels: 

Cut  102  -  20"  long  X  12"  wide  pieces 

Trim  the  4"  wide  X  40"  long  pieces  into  34  -  20"  long  X  4"  wide  pieces 


Figure  3.6  -  Prepreg  Cutting  Plan  for  90  Plies 


The  baling  sequence  Listed  in  table  3.)  below  and 
illustrated  in  figure  3.7  was  identical  for  ail  of  the 
laminates  manufactured: 

1  .  Vacuum  bag 

2.  Vent  blanket  breather  (Air  Weave  SS-FR) 

3.  Perforated  mylar  (seal  to  coroprene  edge) 

4.  Mochburg  bleeder 

5.  Perforated  mylar  (seal  to  coroprene  dam) 

6.  TX  1040  porous  teflon  coated  glass  cloth 

7.  Eight  ply  laminate 

8.  TX  1040 
y  .  Mylar 

10.  12  inch  radius  of  curvature  steel  mold 

Table  3.1  -  Bagging  Sequence  (3) 


Vacuum  Bag 


Breathe: 


Perforated  Mylar 
B’eeaer 

Perforated  Myar 


TX  tO-O 


Laminate 


TX  1040 


Mylar 


Tcoi  Steel  Mo 


The  laminates  were  then  cured  in  an  autoclave  according 
to  cycle  B-240-T,  Rev  A  as  shown  in  table  3.2  below : 

1.  Apply  Vacuum,  25  in  Hg  and  85  ps l  minimum 

2.  Heat  air  to  240  F  in  30  :  5  min 

3.  Hold  part  at  240  :  5  F  for  60  min  under  85  ps i  and 
full  vacuum 

4.  Increase  pressure  to  100  psi  and  vent  vacuum 

5.  Heat  air  to  350  F  in  30  ;  5 

6.  Hold  part  at  350  !  5  F  and  100  ps:  for  120 
minutes 

7.  Cool  part  below  150  F  in  no  less  than  120  min 

8.  When  part  is  below  150  F,  vent  pressure 

Table  3.2  -  Autoclave  curing  cycle  used  [3] 

After  curing,  the  panels  were  removed  from  the  bagging 
and  visually  inspected.  At  this  time,  it  was  noted  that  one 
mold  pr  uluced  a  panel  which  had  a  glossier  sheen  than  the 
other  panels.  A  resin  analysis  showed  that  these  panels  had 
a  higher  than  normal  void  content,  with  a  range  of  2.0  to 
4.5  perc.-nt  voids.  Normally  these  panels  would  have  been 
discarded,  but  it  was  decided  to  test  them  to  determine  the 
effei  ts  of  this  manufacturing  defect  on  panel  strength. 


After  subsequent  autoclave  runs,  it  was  discovered  that 
one  of  the  autoclave’s  vacuum  lines  was  partially  blocked 
and  was  not  producing  a  sufficient  vacuum,  thus  resulting  in 
the  higher  void  contents  in  the  affected  panels. 

After  removal  from  the  molds,  the  panels  were  C-scanned 
using  a  hand  scanner  to  insure  that  no  large  voids  or 
delaminat  ions  were  present.  During  the  L'-scanning  process, 
the  inserted  de 1 am i na t i on s  were  marked  so  that  the  panels 
:ould  be  accurately  trimmed.  The  panels  were  then  checked 
at  13  locations  to  determine  the  average  panel  thickness  and 
the  average  ply  thickness.  Thickness  variation  within 
individual  panels  was  small  (on  the  order  of  0.3%);  however, 
the  average  panel  thickness  variations  from  panel  to  panel 
imongst  the  group  of  panels  tested  was  6.8%. 

For  a  given  panel  size  and  a  material  with  constant 
properties,  the  only  variability  that  will  affect  the 
'Tit.ieal  buckling  load  is  the  panel  thickness ,  h. 
Therefore,  for  panels  with  an  aspect  ratio  of  one,  panel 
buckling  strength  should  vary  by  the  ratio  [  1): 

1'2  -  P  (  h,?/h  (  ) n  (8.1) 

1  r  :  i  a  ’  [.  I  at  n  is  equal  to  three;  but  for  the  curved 
pan  -  .  s  used  in  •'Tils  stud>,  a  value  of  n-8  gave  collapse 
toads  that  diverged  from  STAGSC-  1  pred  i  ct.  ions  . 


However,  from  equation  (3.1),  the  buckling  load  for 
panels  is  seen  to  be  sensitive  to  the  thickness  of  the 
panel.  This  variation  can  cause  scattering  of  the  test  data 
and  experimental  strengths  that  are  other  than  the 
predicted  values.  Equation  (3.1)  was  applied  to  Horban ' s 
STAGSC-1  buckling  loads  for  curved  panels,  and  an  empirical 
regression  relation  was  derived  for  panels  with  a  13  inch 
radius  of  curvature  (see  fig  3.8): 

P2  =  PY  ( h  2 / h 1  )2‘°6  (3.2) 

This  relation  proved  very  useful  in  this  thesis  since  a 
typical  STAGSC-1  linear  bifurcation  analysis  used  from  120 
to  150  cp  seconds  of  computer  time.  Using  equation  (3.2), 
the  largest  calculated  error  was  0.86  percent.  Therefore,  a 
STAGSC-1  run  was  not  necessary  for  every  panel  thickness  in 
the  test  group,  and  the  number  of  STAGS  computer  runs  was 
substantially  reduced. 

Next,  two  test  panels  were  cut  from  each  mold  panel 
using  a  curved  steel  fixture  to  hold  the  panel  and  a  radial 
arm  saw  with  a  water  lubricated,  seven  inch,  diamond  tipped 
blade.  The  panels  were  trimmed  to  a  13  inch  height  by  a 
inch  arc  length  to  allow  for  mounting  in  the  test 
fixture.  The  unsupported  dimension  of  the  panel  when  in  the 
fixture  was  12  inches  by  12  inches. 
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Compression  Load  (pounds) 


Figure  3.8  -  STAGSC-1  Regression  Curve  for  12"  Curved  Panels 
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Laminate  thickness 


Strain  gauges  (  M&M  CEA- 2  5  0- 3  50UW- 3  50  )  were  placed  on 
each  panel  to  determine  the  uniformity  of  load  introduction 
to  the  panel  and  to  determine  the  load  at  which  strain 
reversal  occurred  at  the  delamination  (fig  3.9).  Strain 
gauges  1  and  2  (SGI  and  SG2 )  were  placed  on  the  concave  and 
convex  sides  of  the  panel  respectively.  These  gauges  were 
at  the  center  of  the  panels’  damaged  area  so  that  moments, 
eccentricities  and  strain  differentials  due  to  pane  1  damage 
could  be  determined.  This  panel  damage  is  seen  as  a  strain 
reversal  between  SGI  and  SG2  as  moments  due  to  panel  damage 
cause  bending  in  the  damaged  region  of  the  panel.  Gauges 
SG3  and  SG4  were  placed  at  the  top  edge  of  the  panel  to 
track  the  uniformity  at  which  the  top  edge  displaced  during 
the  loading  process. 

From  previous  work  [6,7,12],  it  was  determined  that  the 
maximum  tolerance  in  panel  height  should  be  held  to  0.01 
inches  to  ensure  an  even  load  distribution  within  the 
panel.  When  this  tolerance  is  exceeded,  non-symme trie 
buckling  may  result.  Horban  found  that  the  machine  shop 
that  cut  his  panels  was  able  to  trim  his  panels  to  within  3 
mils  of  the  prescribed  length.  Variations  of  the  arc  length 
dimension  weren’t  as  critical  since  a  simple  knife  edge, 
vertical  support  was  used,  and  since  no  transverse  load  was 
appl i ed . 

Scrap  pieces  left  over  after  trimming  were  tested  for 
void  content,  resin  content,  volatile  content,  and  density. 


content 


Void  content  is  an  excellent  indicator  of  panel 
uniformity  (see  table  A.l).  Therefore.,  panels  were  tested 
which  had  high  void  contents  to  determine  the  effects  of 
this  defect  on  panel  strength.  Croop  [3]  wrote  a  report 
detailing  these  deviations  in  panels  that  Horban  |7|  and 
Siefert  [12]  had  used  for  their  research.  Based  on  the  data 
collected  by  Croop,  it  was  assumed  that  panels  with  void 
contents  much  above  0.8%  had  suffered  degradation  during 
panel  manufacture,  and  their  material  properties  were 
suspect . 
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Prepreg  -  Hercules  12  inch  tape,  AS4/3501-6 

Release  Agent  -  RAM  mold  release  225,  a  tooling  release 

manufactured  by  RAM  chemicals 

Teflon  Inserts  -  DuPont  Teflon  FEP  flourocarbon  film, 

0.5  mil  300A  film 

Strain  Guages  -  Micro  Measurements  CEA-250LW-350 

uniaxial  strain  guages 

Load  Cell  -  Intarface  model  1220-BF,  25k  pound 

Control  Panel  -  MTS  model  413.05 

MTS  406  controller 
MTS  436  control  unit 
Hewlett  Packard  3450A 

Amplifiers  -  NEF  126 

Signal  Conditioners  -  PPM  model  SG-15-2 


Table  3.3  -  Experimental  Equipment  List 


Experimental  Setup 


The  test  fixture  used  in  this  thesis  is  a  modification 
of  the  original  General  Dynamics  1974  design  that,  was  us**d 
by  Wilkens  [20]  in  his  early  work  with  curved  composite 
panels.  Since  curved  fixtures  of  this  type  are  rare, 
experimental  results  should  be  particularly  valuable  to 
designers  working  with  curved  panel  geometries. 

The  panel  test  device  consists  of  a  25,000  pound 
capacity  MTS  compression  machine  with  an  Intarface  model 
1220-BF  load  cell  (see  material  list  -  table  3.3).  The 
panel  is  mounted  in  a  fixture  (fig  3.11)  which  clamps  the 
top  and  bottom  edges  of  the  panel,  and  applies  a  knife  edge 
support  with  negligible  restraining  moments  to  the  side 
supports.  These  boundary  conditions  can  be  thought  of  as 
representative  of  the  edge  restraints  caused  by  a  continuous 
panel  restrained  against  rotation  at  one  end  by  a  major 
structural  support  while  being  restrained  by  panel 
stiffeners  in  the  transverse  direction. 

The  clamped  conditions  on  the  top  and  bottom  edges  are 
achieved  by  using  the  machined  steel  plates  shown  in  figure 
3.12.  and  a  series  of  six  steel  chucks  which  are  machined 


so  that 

they 

have 

a  12 

inch  radius  of 

curvature 

at  their 

interface 

with 

the 

curved 

pane  1  . 

The 

chucks  are 

1/2  inch 

high  and 

are 

held 

firmly 

against 

the 

panel  using 

two  set. 

screws  per  chuck. 
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The  vertical  supports  (tig  3.13)  are  designed  to  apply 
a  knife  edge  support  which  will  allow  displacement  in  the  x 
and  y  directions  and  rotations  in  the  x , ^  direction  while 
restraining  displacements  in  the  z  direction.  The  rubber 
"o”  rings  cause  the  knife  edge  supports  to  place  a  constant 
normal  force  on  each  panel.  A  piece  of  teflon  tape  was 
attached  to  both  sides  of  each  panel  tested  to  further  aid 
in  the  panel’s  freedom  to  move  in  the  x  and  y  directions. 

The  fixture  gives  the  following  edge  supports: 


Top  edge:  v=w=w,  =0  u^free 

Bottom  edge:  u=v=w=w,  =0 

x 


Vertical  edges: 


w  =  w,  =0  u=v=w,  =free 
x  y 


The  fixture’s  vertical  supports  allow  1/4  inch  of  the 
panel  to  protrude  above  the  top  of  the  supports,  to  allow 
for  loading  head  displacement.  An  aluminum  flange  located 
at  the  center  of  the  base  plate  provides  a  mount  for  an  LVDT 
(Linear  Variable  Differential  Transducer)  which  measures  the 
vertical  displacement.  And,  since  the  bottom  test  fixture 
is  not  permanently  attached  to  the  compression  machine,  this 
flange  also  provides  a  mount  to  secure  the  bottom  test  head 
in  a  uniform  position  for  each  test  (fig  3.14).  An  array  of 
15  LVDT ’ s  was  used  to  measure  relative  movement  of  the  panel 


in  '■fie  z  d  i  rec  t  i  or.  dur  i  ng  pane  1  bu>  -k  1  i  ng  (fig  3.15)). 


Test  Procedure 


After  the  panel  was  placed  in  the  test  fixture,  a  30U 
pound  seating  load  was  applied  and  all  of  the  supports  were 
tightened.  The  purpose  of  this  seating  load  was  to  insure 
that  the  panel’s  top  and  bottom  edges  are  clamped  paraLlel 
to  the  heads  and  to  reduce  the  chance  of  introducing  a 
moment  load  at  the  supports.  Tests  with  previous  panels 
with  no  seating  load  indicated  a  large  displacement  was 
induced  before  the  panel  began  taking  load  (fig  3.16).  The 
displacement  was  actually  the  seating  of  the  panel  in  the 
clamped  supports.  The  seating  load  was  then  removed,  and 
all  of  the  channels  were  set  to  zero. 

The  panels  were  then  compressively  loaded  with  a 
constant  0.05  inch  per  minute  displacement  with  a  load  cell 
measuring  the  applied  load.  The  displacement  was  introduced 
at  the  top  edge  of  the  panel  and  is  measured  using  a  LVDT 
which  measures  the  relative  movement  of  the  top  head  with 
respect  to  the  bottom  head  of  the  compression  machine.  A 
matrix  of  15  LVDT ’ s  was  positioned  normal  to  the  convex  side 
of  the  panel  to  measure  deflections  in  the  z  direction 
during  panel  loading.  These  LVDT ’ s  were  positioned  at  x=3, 
6,  and  9  inches  with  the  LVDT ' s  equally  spaced  along  the 
rows  . 

Data  was  recorded  three  times  per  second  for  the  16 
LVDT’s,  the  load  cell,  and  the  four  strain  gauges  with  each 


instrument  assigned  to  its  own  channel.  The  data  was  saved 
on  a  VAX  11/780  computer  and  later  to  magnetic  tape  backup 
so  that  experimental  graphs  could  be  plotted.  In  addition, 
a  CRT  was  used  which  allowed  the  channels  to  be  monitored 
during  loading  so  that  the  progression  of  panel  loading 
could  be  tracked  during  the  test. 

When  global  buckling  of  the  panel  occurred,  the  load 
was  put  on  hold  and  the  buckled  regions  were  outlined  with  a 
silver  pencil  and  photographed.  The  load  was  then  removed 
by  releasing  the  hydraulic  pressure. 

Ten  of  the  panels  tested  had  a  3000  pound  seating  load 
applied  prior  to  testing.  This  large  seating  load  was  used 
to  determine  if  there  was  a  difference  in  the  strain 
reversal  at  the  delamination  as  a  result  of  ply  separation 
brought  about  by  the  high  initial  load.  This  seating  load 
had  the  effect  of  smoothing  the  introduction  of  the  test 
load  since  there  was  very  little  initial  displacement  of  the 


panel  in 

the 

supports 

at 

initial 

test  1 oads  (fig 

3.17). 

However , 

this 

seating 

load 

had 

no 

apparent  effect 

on  the 

1 uad  a  t 

which 

snapping 

occurred 

and 

didn ’ t  seem  to 

effect 

the  shape  of  the  strain  separation  curve.  By  observing  this 
curve,  it  can  also  be  seen  that  the  assumption  of  linear 
material  properties  arid  linear  bifurcation  is  an  accurate 


a ppm\ i ma t  ion. 
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IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


Panel  Identification 

A  panel  identification  nomenclature  similar  to  Horban's 
[7]  was  adopted  for  this  thesis  and  is  shown  in  table  4.1 
below.  Tabulations  of  experimental  results  can  be  seen  in 
tables  5.3  and  A.l. 


Q2 -12-1** 

i 

—  One  asterisk  indicates  3000  lb  seating  load 

-  Two  asterisks  indicate  2  pieces  of  teflon 

used  to  create  delamination 

- Test  number,  panel  number  1,  2,  3  etc. 

- Delamination  located  between  these  plies 

-  Diameter  of  delamination  in  inches 

-  Quasi-isotropic  panel 

Table  4.1  -  Panel  identification 


Panels  with  varying  levels  of  damage  were 
experimentally  tested  by  applying  a  uniform  rate  of  end 
displacement.  The  panels  tested  for  this  thesis  had  damage 
due  to  faulty  manufacturing  of  the  panels,  surface  damage 
similar  to  that  caused  by  burning  or  chafing,  or  an  internal 
delaminat.ion  of  the  type  that  might  occur  as  a  result  of  a 
low  velocity  impact.  Other  panels  with  no  damage  were 
tested  as  a  control  group.  The  panels  were  loaded  until 
global  panel  buckling  occurred,  and  the  buckling  load, 
buckling  strain  and  the  buckled  shape  of  the  panels  were 
noted.  The  responses  of  the  panels  in  their  damaged  regions 
were  recorded  both  visually  and  through  the  use  of  strain 
gauges  centered  on  either  side  of  the  damaged  region. 


Panels  with  No  Damage  -  QO-OO-x 

Panels  which  had  no  surface  damage  and  no  internal 
delaminations  were  experimentally  tested.  These  panels  all 
had  void  contents  within  an  acceptable  limit  of  one  percent 
and  exhibited  no  discontinuities  when  C-scanned.  Therefore, 
it  was  assumed  that  no  manufacturing  induced  defects  were 
present  .  This  group  of  control  panels  had  strain  gauges  SGI 
and  SG2  positioned  at  the  center  of  the  panel  in  the  same 
manner  as  for  the  damaged  panels.  These  strain  gauges 
indicated  a  uniform  rate  of  strain  at  the  center  of  t.he 
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panel  with  no  difference  in  strain  between  SGI  on  the  first 
ply  and  SG2  on  the  eighth  plv  (fig  4.1).  This  uniform  state 
of  strain  was  interpreted  as  a  state  of  axial  strain 
without  bending  in  a  vertical  strip  of  the  panel  through  the 
panels  center.  This  approximation  is  valid  since  the  strip 
is  flat  in  the  longitudinal  direction  and  has  a  relatively 
flat  radius  of  curvature  and  is  not  restrained  in  the 
transverse  direction.  If  moments  or  eccentric  loading  had 
been  present,  their  effects  would  have  shown  up  as  a 
differential  state  of  strain  between  SGI  and  SG2 .  This 
uniform  strain  was  observed  for  each  of  the  QO-OO-x  panels 
tested. 

The  panels’  global  buckling  values  were  recorded  and 
compared  to  a  STAGSC-1  finite  element  model  prediction.  The 
experimental  values  were  lower  than  STAGSC-1  predictions  by 
a  mean  of  12.9  percent  with  a  standard  deviation  of  2.8  (fig 
4.2).  The  fact  that  predicted  buckling  loads  were  somewhat 
higher  than  experimental  loads  was  expected,  since  for 
cylindrical  shells,  a  very  slight  surface  imperfection  can 


result 

in  a 

large 

decrease 

i n  buck  1 i ng 

loads  . 

Also  , 

the 

STAGSC-  1 

mode  1 

was 

a  1  i  nea  r 

bi furcat ion 

model 

wh  i  ch 

was 

expected 

to 

yield 

higher 

buck  1 i ng 

values 

than 

the 

experimental  tests. 
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Figure  4.2  -  Global  Buckling  for  QO-OO-.v  Panel 


Avg  Diff  -  12.9* 
Std  Dev  -  2.9* 
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Manufacturing  Induced  Defects 


When  the  first  sets  of  panels  manufactured  for  this 
thesis  were  removed  from  the  autoclave,  it  was  noted  that 
one  of  the  four  molds  produced  a  panel  which  had  a  glossy 
sheen  as  stated  previously.  The  panel  appeared  norma L  in 
every  other  respect  and  there  were  no  other  visible  signs  of 
defects.  This  panel  was  C-scanned  along  with  the  rest  of 
its  group,  and  when  no  voids  or  cracks  were  detected,  it  was 


t  r i ramed 

to  it. 

s  test  dimensi 

ons.  The 

pieces 

trimmed  from  the 

pane  1  s 

were 

analyzed  for 

density, 

fiber 

content  ,  resin 

content 

and 

void  content 

(see  ref. 

.  3  )  . 

It  was  found  that 

the  panel  with  the  sheen,  had  a  much  higher  than  normal  void 
content  (on  the  order  of  2  to  4  percent  as  opposed  to  the 
less  than  1  percent  void  content  for  a  "normal”  panel). 

In  the  mean  time,  other  autoclave  runs  had  produced 
panels  with  the  same  glossy  appearance,  all  from  the  same 
mold.  When  a  lab  analysis  of  these  panels  was  carried  out, 
it  was  found  that  they  also  had  a  higher  than  normal  void 
content.  At  this  time  the  autoclave  was  investigated,  and 
it.  was  found  that  one  of  the  sets  of  vacuum  lines  was 
partially  blocked  with  resin,  and  wasn’t  pulling  sufficient 
vacuum  on  the  bagged  samples.  Apparently,  this  resulted  in 
the  incomplete  removal  of  volatiles  during  panel  curing 
causing  inadequate  consolidation  of  the  panels  and  higher 
void  contents.  After  the  autoclave  was  repaired,  ai 1  panels 
had  a  void  content  below  one  percent . 
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Since  the  number  of  panels  allowed  for  this  research 
was  limited,  it  was  decided  to  test  the  panels  with  high 
void  contents  and  observe  their  buckling  characteristics. 
These  panels  were  assumed  to  have  the  same  material 
properties  as  low  void  panels.  They  were  compress  1 ve 1 y 
loaded,  and  it.  was  found  that  they  had  buckling  loads  which 
were  within  a  mean  of  9.6  percent  of  their  STAGSL'-l  values. 

This  result  was  unexpected.  It  was  thought  that  a  high 
void  content  was  an  indicator  of  a  faulty  specimen  and  that 
these  panels  would  have  greater  flexibility  due  to  the  voids 
which  would  result  in  reduced  buckling  loads.  When  percent 
voids  was  plotted  against  compressive  load  at  failure  and 
the  best  fit  line  was  plotted  (fig  4.3),  it  was  found  that 
strength  was  increasing  as  void  content  increased.  This 
phenomenon  was  explained  by  checking  the  average  thickness 
of  the  panels.  Voidy  panels  were  thicker  than  "normal'' 
panels  (fig  4.4),  and  the  increase  in  thickness  was 
stiffening  the  panels  and  thus  compensating  for  the  higher 
void  contents.  As  a  result,  the  panels  with  a  higher 
percentage  of  voids  had  a  higher  average  buckling  strength 
than  the  panels  with  a  normal  void  content.  However,  in 
general ,  high  void  contents  should  be  avoided  in  a  composite 
panel.  As  the  void  content  increases,  the  amount  of  matrix 
bonded  to  each  fiber  decreases  which  results  i  ri  a  reduction 
of  load  transfer  through  shear  stress.  Moisture  absorption 
and  subsequent.  material  degradation  could  also  be  expected 

to  increase. 
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Surface  Danage 


i 

i 

Surface  damage  in  a  composite  panel  may  occur  to  the 
exterior  surface  or  to  the  interior  surface  of  the  panel  as 
the  result  of  burning,  chafing  or  abrasion.  These  tvpes  of 

i 

surface  damages  were  modeled  in  this  thesis  by  the  removal  I 

I 

i 

ot  a  circular  portion  of  an  exterior  0  ply.  The  panels 
were  then  tested  experimentally  to  determine  the  effects  of 
this  damage  on  the  panels'  buckling  strength. 

These  Q2-12,  Q2-78,  Q4-12  and  Q4-78  panels  were 

manufactured  with  a  teflon  disc  placed  between  plies  1-2  and 
between  plies  7-8.  The  de 1  am i nat i ons  caused  by  placing 
inserts  adjacent  to  the  exterior  ply  resulted  in  the 
de 1 am i na t i ons  that  were  clearly  failed  (fig  3.1).  Since  the 
blistered  portion  of  the  external  ply  was  failed  during  the 
manufacturing  process,  it  was  removed  and  SGI  and  SG2  were 
placed  at  the  center  of  the  damaged  region  on  the  exterior 
of  the  90  and  the  45  plies.  These  strain  gauges  were  used 
to  detect  differential  strain  in  the  damaged  region  which 
would  indicate  the  presence  of  eccentric  loading  and  moments 
due  to  the  damage.  It  was  thought  that  the  additional 
moments  caused  by  this  eccentricity  would  trigger  buckling 
at  lower  values  than  those  observed  for  nori-de  1  am  l  na  t  eo 
panels.  But,  the  experimental  buckling  loads  differed  from 
H."  STAGSC-1  predictions  by  a  mean  of  only  tt  .  K  per  ••n'  with 
a  standard  deviation  of  5.8  percent  ifig  1.5..  !'he  vis  a 


observable  difference  be  tween  panels  with  surface  damage  < .  n 
the  ronvex  side  as  opposed  to  panels  with  damage  to  the 
concave  surface.  The  buckling  values  for  these  panels  were 
actually  higher  buckling  loads  than  were  recorded  for  the 
non -de 1  am i na ted  panels  and  are  very  close  to  the  values 
recorded  for  panels  with  high  void  content.  This  fact  may 
be  explained  by  noting  that  very  thin  eight  ply  panels  were 
used  for  this  thesis.  All  surfaces  are  very  close  Lo  the 
neutral  axis;  and  therefore,  bending  moments,  which  change 
by  the  distance  from  the  neutral  axis  cubed,  are  small.  As 
a  result,  eccentricities  caused  by  this  surface  damage  are 
minuscule;  however,  for  thicker  panels  with  deeper  surface 
damage,  these  moments  could  increase  and  could  cause  reduced 
panel  buckling  loads.  The  total  volume  of  panel  that  is 
removed  as  a  result  of  this  surface  damage  is  also  a  very 
small  percent  of  the  total  panel  volume.  For  a  four  inch 
iiamet er  surface  damage  in  the  outer  ply  of  the  panel,  only 
1.1  percent.  of  the  panel’s  volume  is  removed.  As  a  result  , 
the  panels  are  seen  to  respond  locally  to  surface  damage, 
while  little  effect  is  seen  on  global  panel  buckling. 

The  local  effects  of  this  damage  can  be  seen  by 
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Hoover,  when  the  Q4-1H  and  w 4-78  panels  were  tested, 
an  opposite  state  of  relative  strain  existed.  For  these 
panels,  the  strain  gauge  on  the  damaged  surface  was  placed 
in  tension  relative  to  the  gauge  on  the  undamaged  surface. 
Also  noted  was  a  dramatic  increase  in  the  relative  strain 


differential  at  the  center  of  the  damaged  region  at  global 
panel  buckling  (fig  4.6). 

This  strain  reversal  could  tie  explained  if  the  damaged 
region  were  to  be  thought  of  as  a  plate.  As  compressive 
load  is  applied,  a  load  is  reached  which  triggers 
instability  of  the  damaged  region.  The  damaged  region  is 
restrained  from  buckling  by  the  surrounding  panel  which  acts 
as  a  clamped  support,  but  it.  eventual  ly  readies  a  lead  at 
which  1  oca  1  buckling  occurs.  For  larger  damaged  regions, 
such  is  the  four-  inch  diameter  abrasion,  a  second  mode 
buck! i ng  miy  occur  causing  a  reversal  of  r*  i a t : \ e  strain  on 
a  lire  draw  n  'hru  the  damaged  region  and  thus  tension  at  t  tie 


e  n  t  i 

.*  r 

(  %  i 

t  h  e 

damaged  area. 

More  e 

■  x  p>e  r  i  me  tit  al  or  1 

i  u  i 

e  i  e  m <■ 

nt 

W  <  t 

rk 

o u ! d  have  >o  be 

doin'  t 

•  determine  t  he 

r  i 

s  *  a  t  > 

*s 

i  n 

t  t.  e 

damaged  region 

a  nd  t  <  t 

•‘heck  t  h"  valid] 

i  r  > 

this 

by ; 

lot  h 

•  •  SIS. 

6  H 


a  6000 


4000 


-  ”  "  A^  * 

A  A  4  A 

B 

□ 

A 

-  STAGSC— 1  Prediction 

*  Q2-12  Panels  Avq  Diff 

0  Q2-78  Panels  Std  Dev 

9. a* 

5.6* 

Laminate  Thickness  (mils) 


'■\d  vs  Laminate  Thickness  -  Surface  Damaet 


>0  001 


Internal  Delamination 

Impact  damage  to  a  composite  shell  is  usually  in  the 
form  of  an  internal  delamination.  This  type  of  damage  was 
modeled  using  teflon  disc  inserts  to  prevent  bonding  between 
adjacent  plies.  In  previous  work  [7],  it  had  been  found 
that  some  insert  materials  didn’t  create  total  debonding. 
Therefore,  to  determine  if  100  percent  delamination  at  the 
inserts  was  occurring,  a  panel  was  tested  by  C-scanning  and 
by  using  stereo  x-ray.  The  C-scan  indicated  that  a 
discontinuity  existed  at  the  approximate  boundaries  of  the 
delamination.  The  stereoscopic  x-ray  confirmed  that  100 
percent  delamination  was  being  achieved. 

The  delaminated  area  of  one  of  the  panels  is  shown  in 
figure  4.7.  This  stereo-scopic  x-ray  was  produced  by 
introducing  a  radioactive  isotope  penetrant,  tetrabro- 
methane,  into  the  panel  through  a  3/8  inch  hole.  After 
approximately  half  an  hour,  the  penetrant  had  worked  its  way 
'  la  capillary  action  into  all  accessible  cavities.  These 
areas  show  up  as  light  areas  on  the  photograph.  By  looking 
closely,  an  observer  can  easily  see  the  outline  of  the 
circular  delamination,  and  can  determine  that  complete 
delamination  has  occurred.  This  test  was  only  done  on  one 
panel  since  extensive  work  had  been  done  previously  bv 
Horban  [7]  to  determine  the  extent  of  de 1 am i na t i ons  in 
composite  panels. 


Panels  with  two  inch  delaminations  were  then  tested 
experimentally.  These  Q2-23  and  Q2-67  panels  experienced 
global  panel  buckling  at  a  load  which  was  a  mean  of  13.2 
percent  lower  than  STAGSC- 1  predictions  with  a  standard 
deviation  of  5.5  percent  (fig  4.8).  Similar  panels  tested 
by  Horban  had  average  buckling  loads  that  were  19.6  percent 
lower  than  the  predicted  loads.  These  experimental  buckling 
loads  for  delaminated  panels  were  only  slightly  lower  than 
the  experimental  buckling  loads  for  undamaged  panels. 
However,  initial  sublaminate  instability  was  detected  at  a 
much  lower  load.  When  these  panels  were  compress ively 
loaded,  the  delaminated  region  (the  sublaminate)  would  reach 
its  local  critical  load  and  buckle  or  ''snap".  This  snapping 
effect  can  be  detected  visually  as  a  blistered  region  at  the 
delamination  (fig  4.9  and  4.10)  or  it  can  be  detected 
experimentally  through  the  use  of  SGI  and  SG2  centered  over 
the  delamination. 

This  local  sublaminate  snapping  causes  a  release  of 
strain  energy  and  a  lessening  of  the  compression  on  the 
strain  gauge  attached  to  the  sublaminate.  A  strain  reversal 
effect  results,  an  ideal  depiction  of  which  can  be  seen  in 
figure  4.11.  In  this  figure,  the  strain  reversal  is  seen  to 
occur  at  approximately  1900  pounds.  At  local  sublaminate 
buckling,  SGI  on  the  base  laminate  is  in  compression,  while 
SG2  on  the  sublaminate  is  seen  to  have  its  compressive 
strain  relieved.  SG2  goes  into  tension  as  the  delaminated 
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Figure  4.8  -  Experimental  Load  vs  STAGSC-1  Load 
for  Delaminated  Panels 


plies 


snap  out. 


Another  effect  that  can  be  seen  by 


observing  figure  4.11  is  the  lessening  of  compression  on  the 
base  laminate  at  approx  1 mate  1 y  3700  pounds.  it  is  thought 
that  this  corresponds  to  an  initial  panel  deformation  just 
prior  to  global  panel  buckling.  An  attempt  was  made  in 
this  thesis  to  predict  this  sublaminate  instability  through 
the  use  of  the  linear  superposition  of  simple  finite  element 
models  and  a  plane  strain  adjustment.  This  method  will  be 
presented  in  chapter  five. 

During  the  experimental  testing  of  the  panels,  strain 


c 

was 

o  n  1  y 

obse rved 

in  the  axial 

direction.  Also,  an 

axial 

load 

was 

appl ied 

to  the  panel 

in 

the  plane  of  the 

panel  . 

Th  i  s 

led 

to  the 

observation  that 

a  state  of  plane 

strain 

could  be  approximated  for  the  panel  in  a  strip  in  the  region 
of  the  de lam i nat i on . 

The  delaminated  region  of  the  panel  was  closely 
O  observed  after  ex per imen  tal  testing.  This  region  was 

i  observed  to  have  undergone  a  plastic  deformation  as  a  result 

I 

of  panel  buckling.  When  the  panel  was  unloaded  and  removed 
%  from  the  test  fixture,  the  delaminated  area  was  clearly 

[  visible,  and  retained  the  same  general  deformed  shape  it  had 

attained  under  buckling  load  (fig  4.12).  This  indicated 
that  this  region  had  undergone  plastic  deformation  and 
indicated  a  possibility  of  delamination  growth. 

At  this  point,  a  number  of  panels  were  placed  back  into 
the  test  fixture  and  compres s i ve 1 y  loaded  for  a  second 
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time.  The  expected  strain  reversal  ot'  SGI  vs  Si,_  i. 
occurred,  but  at  a  much  lower  value  of  load  than  I <  -  r  '?• 
initial  test.  Figure  4.13  shows  that  local  su  t>  I  tn,  .  r  a  • 
buckling  is  occurring  at  a  load  of  approx  1  riia t  e  1  \  loo  n<  o 
as  opposed  to  the  1900  pound  load  seen  for  the  initial 
of  the  panel  (fig  4.11).  This  lowering  of  t  h«-*  buck  .  t ;  g 

load  for  the  delaminated  plies  was  expected,  since  a  plast  i  ■ 
deformation  of  the  sublaminate  would  reduce  its  stiffness, 
would  substantially  increase  the  load  eccentricity  due 
out  of  plane  deformations,  and  would  cause  the  second  order- 
rotational  terms  to  become  more  dominant  to  local  buckling. 
As  a  result,  the  critical  load  for  local  buckling  would  be 
decreased  for  subsequent  panel  loadings. 

When  panels  were  reloaded  after  their  initial  buckling 
failure,  a  small  decrease  in  global  panel  buckling  load  was 
also  noted.  It  can  be  seen  that  the  panel  in  figure  4.13, 
which  has  undergone  plastic  deformation,  buckles  at  33U0 
pounds,  as  opposed  to  the  4100  pound  load  required  to  buckle 
the  undamaged  panel  shown  in  figure  4.11.  This  decrease  in 
global  panel  strength  would  occur  as  a  result  of  the 
increased  instability  of  the  delaminated  region  which  has 
suffered  permanent  deformation.  It  is  postulated  that 
subsequent  cyclic  panel  loadings  would  cause  delamination 
grow'th  at  the  delamination  crack  tip  and  would  cause  a 
gradual  deterioration  of  the  panel’s  ability  to  carry  load. 
This  was  beyond  the  scope  of  this  thesis,  but  is  an 


Figure  4.12  -  Plastic  Deformation  of  Sublaminate 
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important  problem  sini'e 


a  panel  in 


an  aircraft  could  t>e 


expected  to  undergo  cyclic  loadings  prior  to  detection  of 
delamina t inn  damage. 

Another  phenomenon  that  was  noted  was  the  shape  of  the 
buckled  sublaminate.  Most  of  the  two  inch  de  laminat ions 
would  snap'  out  into  a  single  blistered  region  that  appeared 
to  be  a  low  bubble  on  the  surface  of  the  panel.  The  four 
inch  sublaminates ,  on  the  other  hand,  would  snap  into  a 
pattern  of  "ridges  and  valleys"  that  were  oriented  at 
approximately  a  45  angle  to  the  1-direction  (fig  4.10). 
This  would  indicate  that  shear  stresses  were  developing  at 
the  sublaminate  boundary,  and  that  shear,  or  a  combination 
of  shear  and  compression  was  causing  buckling  for  the  larger 
diameter  sublaminates. 

It  is  also  thought  that  the  base  laminate,  which  was  in 
contact  with  the  sublaminate,  was  laterally  bracing  the 
delaminated  plies  and  causing  them  to  buckle  at  higher 
modes.  This  phenomenon  would  tend  to  increase  the  critical 
load  at  which  larger  sublaminates  initially  become  unstable 
relative  to  an  unrestrained  sublaminate. 
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ANALYTIC  TECHNIQUE 
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Introduction 

Composite  laminates  subjected  to  a  low  speed  impact, 
such  as  a  dropped  tool  or  a  manufacturing  load,  often  develop 
an  internal  de laminat ion .  De 1 am i na t  i  ons  in  composite  panels 
cause  a  reduction  of  local  stiffness  of  the  panel  plies  in 
the  region  of  the  delamination.  As  a  result,  when  the  panel 
is  subjected  to  loading,  this  area  is  the  first  portion  of 
the  panel  to  become  unstable  and  buckle.  This  local  buckling 
of  the  delaminated  plies  precedes  the  global  buckling  of  the 
panel  and  occurs  at  a  much  smaller  load. 

The  effects  of  this  local  buckling  will  become  more 
important  to  engineers  and  designers  as  composites  are  used 
more  often  in  modern  high-speed  aircraft.  Local  surface 
defects  can  cause  drastic  disturbances  in  air  flow,  and  a 
weakened  area  of  a  pane  L  may  experience  peeling  and  dramatic 
delamination  growth  as  the  surface  is  exposed  to  increased 
aerodynamic  loading. 

Since  curved  panels  are  3-d i mens i ona 1 ,  and  buckling  is  a 
nonlinear  phenomenon  due  to  the  large  rotations  and  moments, 
the  compressive  load  which  will  cause  curved  panels  to  become 
unstable  is  extremely  hard  to  predict  analytically.  This 


thesis  presents  a  technique  whereby  the  local  buckling  loads 


at  the  lie  1 am i na t 1  on  may  be  predicted  within  approximately  ju% 
ot  experimental  values  by  using  a  2 -d i mens i ona 1  model  with  a 
plane  strain  correction.  This  approach  to  predicting  local 
sublaminate  instability  is  based  upon  the  principle  of  linear 
superposition.  During  the  experimental  phase  of'  this 

researcn,  certain  physical  responses  of  the  panels  led  to 
certain  simplifying  assumptions.  It  was  observed  that  the 
local  instability  and  snapping  of  the  delaminated  plies 
occurred  at  a  much  lower  value  of  load  than  the  global  panel 
buckling  for  the  1 2  inch  curved  panels.  It  was  thought  that 
there  existed  a  value  of  strain  which  when  applied  to  the 

circular  delamination  would  cause  the  sublaminate  to  buckle. 
This  value  of  strain  in  the  two  ply  sublaminate  should  be 

related  to  a  value  of  strain  in  the  eight  ply  curved  panel. 
Some  method  of  relating  these  two  values  of  strain  was 

requ i red . 

Therefore,  a  strip  of  the  panel  in  the  vicinity  of  the 
delaminations  was  considered.  This  strip,  shown  by  the 

dashed  lines  in  figure  5.1,  would  be  flat  in  the  longitudinal 
direction  (which  is  the  direction  of  loading  )  ,  and  s  1  l  gh  ■  1  ■■ 
curved  in  the  transverse  direction.  In  the  transv.- :  > . 

y-direction,  the  strip  is  restrained  against  dpf-ni.v 
the  eight  ply  thickness  of  panel  between  t  h>  • 

fixture's  vertical  support,  while  in  t  ru-  \  c 
the  strip  is  unrestrained  again:-’  :■  ' 
displacement  is  applied  in  »  pc  x -  ■  i  :  r • 
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the  y  direction  and  free  movement  in  the  z  direction. 


Hence , 


a  plane  strain  assumption  is  valid  for  this  region  of  the 
pane  1  . 

For  the  panel  geometry  investigated,  a  load  applied  to 

the  top  edge  of  the  curved  panel  causes  a  displacement  at  the 

boundary  of  the  delaminated  plies.  It  was  found  that  this 

displacement  causes  an  apparent  constant  strain  field  through 

the  curved  panel  in  the  direction  of  load,  and  introduces  an 

effective  strain  to  the  delaminated  plies  via  an  axial 

deflection  of  the  sublaminate  boundary.  The  relation  between 

strain  in  the  curved  panel  and  strain  in  the  delaminated 

plies  was  used  since  strain  is  a  convenient  tool  for 

associating  the  reponse  of  the  curved  panel  with  the  behavior 

of  the  sublaminate  plies.  Since  bending  is  negligible,  the 

delamination  is  primarily  affected  only  by  the  deflection  in 

the  axial  direction  that  is  transmitted  by  a  constant  strain 

in  the  curved  panel.  To  investigate  the  effect  of  this 

deflection  on  the  sublaminate,  a  finite  element  model  of  the 

sublaminate  was  loaded  with  an  edge  deflection.  The  N  load 

x 

and  the  related  axial  strain  in  the  sublaminate  were  computed 

at  local  buckling  .  A  plane  strain  finite  element  model  was 

used  to  relate  this  value  of  N  to  the  eight  ply  curved 

x 

panel.  This  model  represented  the  strip  shown  in  figure  5.1, 
and  has  eight  plies  through  the  z-direction,  and  containes  a 
two  ply  delaminated  region.  This  model  was  loaded  using  the 


strain  which  had  caused  buckling  in  the  circular  finite 
^  element  model  of  the  two  ply  sublaminate. 

The  resulting  value  through  the  cross  section  of  the 
plane  strain  model  was  calculated  at  a  section  of  the  model 
^  one  delamination  radius  away  from  the  delamination  tip.  This 

datum  was  chosen  so  that  the  calculated  force  resultant  would 
be  far  enough  from  the  discontinuity  so  that  model  responses 
0  due  to  the  delamination  tip  could  be  avoided.  This  force 

resultant  in  the  eight  ply  plane  strain  model  could  then  be 
related  to  the  global  panel  load,  and  would  represent  a  force 
resultant  in  the  eight  ply  laminate  corresponding  to  the 
buckling  strain  in  the  (0/-45)  sublaminate. 

In  summary,  the  prediction  of  the  compressive  load  which 
•  will  cause  the  delaminated  region  of  a  composite  panel  to 

snap  was  developed  as  follows,  and  is  outlined  in  Table  5.1: 

1.  A  model  of  the  curved  panel  containing  no 

©  delaminations  was  run  using  the  STAGSC-1  finite  element 

program  for  a  normalized  edge  displacement.  An  N  value  at 
the  top  edge  of  the  panel  was  calculated  for  this 
%  displacement,  and  the  strain  in  the  panel  associated  with 

this  force  was  calculated  at  the  datum  shown  in  figure  5.1. 
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STAGSC-1  CURVED  PANEL  MODEL 

1.  Compute  buckling  load.  P,  lor  non-damaged  panel. 

2  Compute  strain  in  panel  at  global  buckling. 

3.  Use  the  plane  strain  model  shown  below  to  calculate 
a  force  resultant  in  a  strip  of  the  panel  buckling. 
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PUNE  STRAIN  -  NO  DEUMINATION 


t.  Use  the  strain  computed  in  the  curved  panel 
model  above  to  load  the  plain  strain  model 
2.  Compute  the  lorce  resultant  in  this  8  ply  model 
associated  with  the  strain  obtained  in  the 
curved  panel  model  above. 


N 


=  N 


x2 


6 


.  4* 
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STAGSC-1  DISC  MODEL  (0/-45) 


1  This  is  a  model  ol  the  sublammate  created 

by  the  delamination 

2  Compute  the  lorce  resultant  associated 

with  biturcation  tor  this  model 

3  Compute  the  strain  associated  with 

biturcation  lor  this  model 

4  This  strain  will  be  used  to  load  an  8  ply 
plane  strain  model  taken  at  the  delammation 


=  N 


x3 


PLANE  STRAIN  WITH  DEUMINATION 

1  Load  this  model  with  the  strain  calculated 
at  buckling  lor  the  disc  model  above 

2  Calculate  the  lorce  resultant  thru  the  8  ply 

section  associated  with  this  strain 


=  N 


x4 


Predicted  Load  at  Local  Buckling  of  Delaminated  Plies 


pt  -  p 


NX2  Nx3 

N  ,  N  . 

x  1  x  4 


Table  5.1  -  Overview  of  Analytic  Technique 


bifurcation  load  was  determined;  and  the  strain  in  the  disc 
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© 


% 


and  the  N  value  at  plate  buckling  were  determined.  The 
strain  obtained  at  bifurcation  for  this  model  was  then  used 
to  load  the  plane  strain  model  in  step  three. 

3.  A  finite  element  program  [21]  was  used  to  run  a 

plane  strain  analysis  of  an  eight  ply  composite  containing  a 

delamination.  This  model  was  evaluated  using  the  strain 

obtained  at  buckling  for  the  circular  plate  of  step  two. 

This  strain  from  step  two  is  used  to  calculate  an  end 

d i sp lacement  in  the  plane  strain  mode  1  which  ye i Ids  an 

equivalent  strain  at  the  delamination  tip.  For  this  value  of 

strain,  the  value  thru  the  eight  ply  section  of  the  model 

is  calculated.  This  calculation  is  performed  at  a  cross 

section  that  is  one  delamination  radius  away  from  the 

delamination  tip  so  that  stress  concentrations  due  to  the 

delamination  are  avoided,  and  so  that  a  uniform  N  value  can 

x 

be  calculated  which  corresponds  to  an  eight  ply  section  of 
the  curved  panel . 

4.  A  plane  strain  model  was  also  run  for  a  section  of 

panel  that  contained  no  de  lamina  t  i  on .  This  model  was  run 

using  the  strain  in  the  curved  panel  of  step  one;  and  is  used 

to  determine  an  N  value  in  a  plane  strain  model  which  is 

x 

equivalent  to  the  load  applied  to  the  curved  panel.  The  model 


was  constructed  such  that  the  strain  at 

the 

outer  edge 

o  t' 

the 

delaminated  region 

calculated  in  step 

one 

was  used 

as 

the 

input  load  for  the 

plane  strain  model. 

The 

.N  value 

i  n 

an 

88 
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eight  ply  section  of  the  model  located  one  deiamination 


radius  away  from  the  delamination  was  calculated. 

5.  Now,  the  values  have  been  calculated  in  two  plane 
strain  models.  One  force  resultant  corresponds  to  a  uniform 
load  applied  to  the  top  edge  of  a  curved  panel;  and  the  other 
force  resultant  corresponds  to  the  buckling  load  for  a 
(0/-45)  circular  plate.  Using  these  force  resultants,  and 
the  force  resultants  obtained  in  the  two  STAGSC-1  models,  a 
load  is  computed  which  when  applied  to  the  top  edge  of  the 
curved  panel,  will  cause  buckling  of  the  circular 
deiamination.  This  load  is  simply  a  percentage  of  the  load 
applied  to  the  curved  12  inch  panel  (see  table  5.1). 

The  models  used  in  this  analytic  method  and  the  process 
of  calculating  a  local  sublaminate  buckling  load  will  be 
discussed  in  more  detail  in  the  balance  of  this  chapter. 


Curved  Panel  -  STAGSC-1  Model 

One  of  the  standard  geometries  that  STAGSC-1  is  capable 
of  generating  is  a  curved  cylindrical  panel.  For  the  finite 
element  model  used  in  this  analysis,  the  surface  of  the  shell 
was  discritized  into  an  18  by  18  element  mesh  consisting  of 
2/3  inch  square  elements.  Previous  work  done  by  Siefert  I  12 J 
with  this  geometry  had  shown  that  accurate  results  could  tie 
obtained  if  the  element  size  was  between  1/2  and  1  inch.  The 

8y 


element  selected  for  this  analysis  was  the  quadrilateral  410 
element  referred  to  as  the  QUAF  410  (fig  5.2)  in  the  STAGSC-l 
user  manual  [1].  This  element  was  selected  because  it  works 
well  for  linear  analysis  and  for  cylindrical  shells,  and 
because  it  has  12  fewer  degrees  of  freedom  than  the  more 
accurate  QUAF  411  element.  The  QUAF  4  10  is  a  20  degree  of 
freedom  four  noded  element  with  three  translational 
displacements,  two  in-plane  rotations,  and  one  independent 
normal  rotation  at  each  node. 

The  shape  functions  for  the  u  and  v  displacements  within 
the  QUAF  4  10  element  are  a  cubic  polynomial  parallel  to  the 
edges  and  a  linear  polynomial  perpendicular  to  the  edges,  and 
the  bending  shape  functions  are  cubic  polynomials. 

Both  linear  material  properties  and  linear  buckling  were 

used  for  this  model.  The  model  was  loaded  at  its  top  edge 

with  a  constant  end  displacement  in  the  positive  x  direction. 

The  boundary  conditions  used  for  the  model  were  the  same 

constraints  assumed  in  the  experimental  work.  Using  this 

STAGSC-1  model,  the  strains  and  stresses  in  the  panel  could 

be  determined  for  a  given  value  of  load,  P.  Also  computed 

was  the  force  resultant  at  the  top  edge  of  the  panel  ,  N  ,  , 

x  1 

associated  with  this  load. 

The  int.roduct  ion  of  a  load  to  the  top  edge  of  the  curved 
composite  panel  causes  a  strain  m  the  panel  in  the  vicinity 
of  the  delamination  tip.  The  critical  strain  at  this  point 
which  will  cause  the  subiaminate  to  buckle  is  discussed  in 
the  next  section. 


Delaminated  Plies 


STAGSC-1  Model 


A 


f 


The  author  was  interested  in  modeling  the  buckling  of  a 
circular  plate  to  characterize  the  plane  strain  strip  concept 
previously  discussed.  The  most  obvious  condition  was  one  of 
complete  geometric  symmetry,  and  therefore,  the  subsequent 
convergence  study  and  model  was  predicated  on  this  strain 
feature . 

The  delaminations  manufactured  into  the  panels,  as 
stated  previously,  were  two  and  four  inches  in  diameter. 
The  STAGSC-1  model  discussed  in  this  section  was  used  to 
determine  the  force  resultant  and  the  strain  necessary  to 
cause  buckling  of  the  (0/-45)  cicular  delaminated  region  of 
the  panels.  The  force  resultant  at  the  model  boundary,  N  , 

X  o 

and  the  edge  displacement  of  the  boundary  were  determined 
using  a  STAGSC-1  annular  ring  shell  geometry  and  various 
finite  elements  and  mesh  refinements. 

A  convergance  study  was  first  performed  to  determine  an 
optimum  finite  element  mesh  for  this  geometry.  Initially,  a 
mesh  having  12  elements  and  isotropic  material  properties  was 
run  using  the  QUAF  410  element.  The  mesh  was  identified  as  a 
5x4  mesh  since  it  had  five  nodes  on  any  radius  and  four  nodes 
around  any  arc.  The  mesh  was  modeled  as  a  quarter  of  a 
circle  since  symmetric  loading  and  boundary  conditions  could 
be  applied.  A  closed  form  Bessel  function  solution  exists 
tor  a  uniform  axisvmet  ric  compressive  load  on  a  circular 
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plate.  Therefore,  a  uniform  edge  displacement  was  applied 
radially  inward  at  the  outer  boundary,  and  the  value  at 
plate  buckling  was  recorded.  This  model  (fig  5.3)  was  run 
for  both  clamped  and  simple  edge  supports,  and  the  STAGSC-1 
buckling  load  was  compared  to  the  following  Bessel  function 
solutions  for  the  buckling  of  circular  discs  t  4  J  : 


For  clamped  edges  ( u=v=w , x=w , z= f ree  w=w , y=constrained ) : 


14.65  D 

2 

a 


(5.1) 


For  simple  supports  (  u=v=w , y=w , z~ f ree  w  =  w , x=constrained ) : 


where : 


N 


x 


4 . 20  D 

2 

a 


(5.2) 


D  =  „  (o.3) 

a  (  1  ) 

a  =  the  radius  of  the  disc 


It  should  be  pointed  out  that  the  so  called  clamped  edge 
condition  best  characterized  the  actual  delaminated  boundary 
with  respect  to  the  overall  sublaminate  restraint. 

Therefore,  the  boundary  condition  incorporated  in  the 
completed  analysis  is  the  clamped  support. 

The  STAGSC-1  model  was  run  for  increasing  mesh 


refinements  including  5x4,  10x8,  and  20x16  for  both  simple 

and  clamped  edges.  This  study  showed  that  convergence  to  the 


Figure  5.3  -  Meshes  used  for  Ckonvergance 
Study  5x4  and  10x8 


Bessel  function  solutions  was  obtained  with  the  10x8  mesh 
(fig  5.4  and  table  5.2).  During  this  convergence  study, 
initial  models  were  run  with  a  thickness  of  one  inch  to  speed 
the  computations  of  the  closed  form  solution.  These  models 
gave  results  that  diverged  from  the  closed  form  solutions  as 
the  mesh  was  refined  (fig  5.5).  This  is  probably  due  to  the 
fact  that  STAGSC-1  was  developed  using  thin  shell  theory  with 
a  small  linear  rotation  assumption.  Geometries  that  are 
clearly  three  dimensional  (having  a  length  to  thickness 
ratio  less  than  10)  should  be  avoided  when  using  STAGSC-1. 

A  model  used  by  Shivakumar  and  Whitcomb  at  NASA  Langley 
Research  Center  [13]  for  their  research  with  elliptical 
delaminations  was  also  run  (see  fig  5.6).  They  had  obtained 
convergance  with  this  model  for  anisotropic  lay-ups,  so  it 
was  used  in  this  study  to  check  the  accuracy  of  the  10  x  8 
model.  Shivakumar  and  Whitcomb’s  model  uses  the  TRINC  320 
triangular  plate  element  (fig  5.7)  which  has  18  degrees  of 
freedom:  three  translations  and  two  rotations  at  each  of  the 
corner  nodes  plus  edge  rotations  at  midside  nodes.  The 
in-plane  shape  functions  are  linear  for  this  element  and 
bending  is  piecewise  cubic  [10].  This  model  had  738  degrees 
of  freedom  as  opposed  to  the  377  degrees  of  freedom  in  the 
10x8  model  using  the  QUAF  410.  Therefore,  it  was  a  more 
expensive  model  to  run.  Whitcomb’s  model,  using  the  320 
element,  correlated  well  with  the  closed  form  Bessel 
solution,  but  not  as  closely  as  the  models  using  the  QUAF  410 
rectangula  eLeraents. 


Nx  Buckling  Load  (Ib/in  E03)  12  Nx  Buckling  Load  (lb 


Mesh  Refinement  —  Clamped  Edge 

jure  5.4  -  Conversance  for  Circular  Disc  -  h=0.0048 
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Mesh  Refinement  —  Clamped  Edge 

Figure  5.5  -  Divergance  for  Circular  Disc  -  h=1.000" 
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Table  5.2a  -  Circular  Disc  STAGSC- 1  Results 
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Figure  5.6  -  9  x  33  Mesh  with  320  Element 
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Thus,  the  10x8  QlAF  4  10  model  shown  in  ligure  5.7a  was 


•ft 


used  to  obtain  N  at  bucklirug  tor  the  circular  delaminated 

x 

region,  N  Also  determined  using  this  model  was  the 

nominal  f  strain  in  the  sublaminate  at  local  buckling.  This 
strain  was  then  used  in  the  plane  strain  model  discussed  in 
the  next  section  to  determine  an  equivalent  N  value  in  an 
eight  ply  laminate,  N  .  This  force  resultant  is  the  force 
resultant  in  an  eight  ply  plane  strain  model  associated  with 
the  strain  at  buckling  for  a  two  ply,  (0/-45)  circular 
sublaminate . 

Iso  calculated  for  this  STAGSC-1  model  was  the  average 
normal  stress  at  the  sublaminate  boundary  at  local  buckling. 
From  the  STAGSC-1  output  for  a  two  inch  delamination ,  the 
average  N  value  for  an  axisymetric  edge  displacement  of 
0.001  inches  was  48.136  lb/in.  This  value  was  then 
multiplied  by  the  eigenvalue,  0.095,  to  determine  the  average 
force  resultant  at  the  boundary  at  the  onset  of  sublaminate 
instability.  This  force  resultant  was  then  divided  by  the 
cross  section  dimension  of  the  (0/-45)  plies,  2  x  0.0048",  to 
determine  the  average  normal  stress,  476.3  ps  i  ,  at  the 
sublaminate  boundary  at  local  buckling  of  the  model.  This 
stress  will  be  used  in  the  next  section  as  a  comparison  to 
check  the  validity  of  the  plane  strain  model. 
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Plane  Strain  Finite  Element  Model 


0 


0 


The  plane  strain  finite  element  code  used  for  this 
analysis  was  developed  by  Dr.  Sandhu  I  21)  of  the  Air  Force 
Flight  Dynamics  Lab.  The  model  is  used  in  this  analysis  to 
tie  the  strain  at  local  sublaminate  buckling  to  the  strain  in 
a  plane  strain  section  of  the  curved  panel.  The  strain  in 
this  eight  ply  plane  strain  model  could  then  be  used  to 
calculate  a  plane  load  which  would  cause  the  sublaminate  to 
buck le . 

The  problem  was  modeled  as  an  eight  ply  plane  strain 
section  containing  a  delamination  between  the  second  and 
third  plies  (see  fig  5.8).  This  delamination  was  modeled  as 
a  crack  having  no  thickness  in  the  z-direction  by  defining 
two  node  points  at  each  station  along  the  del  am i na 1 1 on 
length.  Also,  the  mesh  is  highly  refined  in  the  vicinity  of 
the  tip  of  the  de lam ination  since  this  is  where  the  greatest 
stress  concentrations  were  expected  to  occur  (fig  5.9). 
Symmetry  is  used  at  the  midlength  of  the  delamination,  with 
this  boundary  fixed  against  rotation  and  displacement  in  the 
x-direction  while  allowing  freedom  in  the  y-direction. 

The  model  was  loaded  with  an  end  displacement 
equivalent  to  the  strain  required  to  buckle  the  (U/-45) 
sub 1 am  mate .  This  displacement  was  applied  one  delamination 
radius  away  from  the  delamination  tip  so  that  a  uniform  cross 
sectional  force  resultant  could  be  achieved. 
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Figure  5.8  -  Plane  Strain  Finite  Element  Mesh  (distorted) 
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F-gure  5  9  Pane  Sfam  Moqe1  at  De^^ma:  cm  I 


I’he  stresses  in  the  vicinity  of  the  del  ami  nation  tip 
were  recorded  and  are  shown  in  figure  5.y.  The  average  ax  1  a i 
stress  in  the  (J/-45  section  of  the  model  was  calculated  from 


these 

stresses  to 

be  4  7  2.1  ps  i  . 

This  value  is 

for  a 

Pi 

strain 

calculat ion 

associated  with 

the  strain  at 

buck  1 

i  n  g 

the  ( 0 / —  4  5  )  STAGSC-1  circular  disc  model.  The  STAGSC- 1 
circular  disc  model  discussed  previously  yielded  an  average 
stress  value  of  476.3  ps l  for  this  same  value  of  strain.  Ihe 
difference  between  these  two  stresses  is  only  il  .  HH  percent  . 
Therefore,  it  appears  that  by  using  a  symmetrical  loading 
pattern,  a  load  distribution  over  the  upper  two  plies  is 
equivalent  to  a  plane  strain  phenomenon.  Also  of  note  in 

this  plane  strain  model  is  the  magnitude  of  st  resses  in  t  ne 
^-direction.  In  the  vicinity  of  the  de 1  am i na !  l  <  >n  tip  these 
stresses  varied  from  U.6H  ps  i  to  1.98  ps  l  .  Phis  level  ,,f 
stress  is  inconsequential,  and  it  can  be  clearly  seen  t  ha'  up 
to  i he  point  of  sub  laminate  buckling,  there  is  no  danger  , 
mode  1  crack  growth  due  to  these  stresses  t  Ihis  was  expo  t  oi 
since  the  loading  is  assumed  i  n  the  p.ane  ot  the 
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Ht  sub  laminate  instability  t  or  a  t  wo  pl  y  (0/-45) 
de 1  am i na  t ion . 

Now,  the  strain  in  the  curved  i  z  incn  panel  resulting 
from  a  panel  load  has  been  related  to  a  force  resultant  in  a 
plane  strain  model;  and  now  the  strain  in  a  circular  disc  at 
local  buckling  has  been  related  to  a  force  resultant  in  a 
plane  strain  model.  It  is  now  possible  to  calculate  the  load 
which  when  applied  to  the  top  edge  of  the  1 2  inch  curved 
panel  will  cause  local  instability  of  a  delaminated  region  of 
the  panel.  This  calculation  is  explained  in  the  next 


Plane  Strain  Model 


No  Delamination 


I 


A  plane  strain  analysis  was  performed  tor  the  laminate 
to  currelate  the  force  resultant  and  strain  in  a  plane  strain 
model  with  the  force  resultant  and  strain  in  a  plane  strain 
strip  of  the  curved  composite  panel. 

the  average  value  distributed  through  the  model  was 
calculated  for  the  strain  (obtained  from  the  curved  panel 
model)  that  was  associated  with  a  displacement  at  the  top 
edge  of  the  curved  12  inch  panel. 

First  the  Q  'values  were  calculated  for  the  material: 
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1  1 


Q 


2  2 


1  -  “l2l'21 


1  -  "l2"21 


(5.4) 


(5.5) 


Q66 


1  2 


(5.6) 


Then  the  values  were  computed  in  the  x-direction , 

where  ^  is  the  angle  between  the  fiber  direction  and  the 
panel’s  x-direction. 


y  j  j  =  Q;  jCos^°  +  2<Q12  +  2y6g  )  sin^cos2^  +  (J^^sin^ 


The  stress  can  be  calculated: 


(o.7) 


IQi;l  ( r  1 


(5.8) 
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The  strain,  c  ,  is  taken  from  the  STAGSL -  1  program.  It 
was  observed,  as  previously  mentioned,  that  the  panel  yields 
a  constant  f  strain.  One  could  quickly  compute  this  strain 
by  using  the  displacement  at  the  panel's  top  edge,  and 
dividing  this  value  by  the  overall  length  of  the  panel,  thus 

yielding  a  of:  fT.x=  ®11  f  (5.9) 

The  stress  is  obtained  for  each  of  the  plies,  and  a 
stress  resultant,  ls  calculated.  This  force  resultant  is 
then  the  value  of  in  the  eight  ply  plane  strain  strip 
associated  with  a  deflection  induced  through  loading  of  the 
curved  pane  1 . 


Analytic  Prediction 


Based  on  the  finite  element  models  discussed  in  the 
previous  sections,  and  using  linear  superposition,  the 
following  equation  was  developed  to  predict  local  buckling  of 
the  sublaminate: 


P  =  p  *  a  *  ’3 


(5.10) 


Where  a  is  the  ratio  of  the  force  resultant  at  buckling 
in  the  STAGSC-1  circular  disc  with  the  force  resultant  in  the 
curved  panel  associated  with  a  load  P  applied  to  the  panel. 
And  ft  is  the  ratio  of  force  resultants  in  the  plane  strain 
models  . 


°  =  Nx3 

/ 

N  . 
x  1 

(  5  . 

,11) 

*  =  Nx2 

/ 

N  A 

x4 

(  5  . 

.  12  ) 

Pj.  =  Predicted  load  in  pounds  at  which  local  sublaminate 
snapping  will  occur. 

P  =  Load  in  pounds  at  top  edge  of  the  curved  panel 

N  ,  =  N  value  at  buckling  for  Curved  Panel 
x  1  x 

N^2  s  Nx  value  in  plane  strain  mode  1  with  No  Delamination 

This  value  is  calculated  for  a  strain  equivalent  to 

the  strain  in  the  curved  panel  due  to  a  load,  P. 

N  _  s  N  value  at  buckling  for  Circular  Disc 
x3  x 

Calculated  using  STAGSC-1  model  of  sublaminate. 

N  ^  =  N'x  value  for  plane  strain  model  containing  a 

delamination.  This  force  resultant  is  calculated 
using  the  strain  at  buckling  for  the  circular  disc 
mode  1 . 


Substituting  equations  (5.11)  and  (o.l2)  into  (5.10), 


the  predicted  load  applied  to  the  top  of  the  curved  12  inch 
panel  which  would  cause  local  buckling  in  the  delaminated 
plies  is : 


N  N 

P  x3  x2 

N  N 

x  1  x4 


( 5 .  13) 


As  the  average  thickness  of  the  laminate  plies 

increases,  the  value  at  which  local  sublaminate  snapping  will 
occur  also  increases.  This  is  seen  by  noting  that  N  (the 
value  for  the  12  inch  curved  panel)  increases  by  panel 
thickness,  h,  raised  to  the  2.06  power  (approximately  by  h^ ) . 
N  ^  the  circular  disc  increases  by  approximately  the  cube 

of  the  panel  thickness;  and  N  _  and  N  .  both  increase 

X  ^  X4 

linearly  with  thickness.  Therefore,  the  snapping  of  the 
sublaminate  was  expected  to  increase  linearly  with  panel 
thickness.  This  linear  variation  of  predicted  snapping  load 
was  plotted  as  a  linear  prediction  of  snapping  load  versus 
panel  thickness  in  figure  5.10.  Using  this  relation,  the 
load  causing  local  instability  of  a  sublaminate  could  be 
determined  based  on  panel  thickness.  This  relation  appeared 
to  be  valid  for  the  range  of  panel  thicknesses  used  in  this 
research . 

The  relations  between  analytic  predictions  and 
experimental  snapping  values  are  shown  in  figure  5.10  and  in 
table  5.3.  The  data  points  shown  in  figure  5.10  are  in 

pairs.  Each  pair  represents  two  panels  from  the  same  mold 


mm 


with  the  same  panel  thickness.  These  data  pairs  are  shown 
connected  by  a  line  with  a  bar  at  midheight  representing  the 


* 
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# 
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average  snapping  load  for  that  pair.  The  snapping  loads  for 
the  two  inch  de  lam  inat  ions  fell  below  the  analytic 
predictions.  This  was  expected  since  linear  models  were  used 
in  the  analysis,  which  in  general  increases  the  load  at  which 
predicted  failure  will  occur.  For  these  panels,  the 
experimental  snapping  loads  differed  from  the  analytic 
predictions  by  a  mean  of  31.3  percent  with  a  standard 
deviation  of  21.4  percent. 

Only  three  data  points  were  obtained  for  the  snapping  of 
the  four  inch  delaminations.  This  is  a  very  sparse  data  set; 
however,  it  was  noted  that  these  panels  tended  to  snap  at 
higher  load  values  than  were  predicted.  A  possible 
explanation  for  this  is  that  the  four  inch  sublaminates  were 
being  restrained  by  the  base  laminate,  and  were  therefore 
buckling  at  higher  modes.  These  panels  differed  from  the 
analytic  predictions  by  a  mean  of  14.6  with  a  standard 
deviation  of  5.2. 

This  analytic  technique  is  seen  to  give  predictions  of 
local  sublaminate  buckling  within  approximately  30  percent  of 
experimental  values.  These  predictions  are  good  engineering 
approximations  obtained  using  a  linear  two  dimensional  model 
to  predict  a  nonlinear,  three  dimensional  local  instability. 
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jure  5.10  -  Strain  Reversal  Load  -  2"  Delaminations 


Figure  5.10  -  Strain  Reversal  Load  -  4"  Delaminations 
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ANALYSIS 


X  DIFF 


Q2-23-1 

40 .69 

1  100 

245  1 

5  5.1 

Q2-23-2 

40.69 

2300 

2451 

6  .  2 

Q2-23-3* 

40 . 80 

1600 

2464 

3  5.1 

Q2-23-4* 

40.80 

2500 

2464 

-1.5 

Q2-23-5* 

41.90 

1000 

2599 

6  1.5 

Q2 - 23-6  *  * 

41.62 

1500 

256  4 

4  1.5 

Q2-23-7  *  * 

42.00 

2200 

261  1 

15.7 

Q2-67-1 

40 . 20 

1900 

2392 

20 . 6 

Q2-67-2 

40.20 

1800 

2392 

24 . 8 

Q2-67-3* 

41  .  10 

2000 

2500 

20 . 0 

Q2-67-4* 

41.10 

1200 

2500 

52 . 0 

Q2-67-5** 

42 . 00 

1000 

261  1 

61.7 

Q2-67-6** 

41.62 

2200 

2564 

14.2 

Q4-23-1 

39.30 

1000 

1168 

14.4 

Q4-23-2 

39 . 30 

1400 

1168 

19.9 

Q4-67-1 

39 . 60 

r>or>r> 

1  187 

9  O  ■>  ■-> 

Q4-67-2 

39 . 60 

1300 

1187 

9 . 5 

NOTES :  * 

-  3000  lb  seating 

load 

tt  -  2  pieces  of  teflon,  3000  lb  seating  load 


Table  5.3  -  Experimental  Delamination  Snapping  Loads 
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VI.  CONCLUSIONS 


Based  on  the  experimental  and  analytic  results,  the 
following  conclusions  can  he  made  for  g r a ph 1 t e / e pox y 
cylindrical  composite  panels  who  h  have  been  subjected  to 
axially  compressive  loads: 

1  .  In  general,  the  linear  bifurcation  predictions  obtained 
using  the  STAGSC- 1  finite  element  code  were  in  close 
agreement  with  experimental  results.  For  undamaged  panels, 
the  experimental  global  buckling  load  was  an  average  of  12.9 
percent  lower  than  predicted  values  with  a  standard 
deviation  of  2.8  percent. 

2.  Surface  damage  to  composite  panels  tends  to  decrease  the 


€* 

panel 

s  t l ffness 

resulting  in  a 

dec  rease 

in  critical 

buck  l  ing 

1  oads  . 

As  the 

surface  area 

that  is 

damaged  i s 

en 1 arged , 

g 1 oba 1 

buck  ling 

loads 

will 

in  general 

dec  rease . 

However , 

% 

this 

study  found  that 

for  2  inch  diamet 

er  surface 

damage  to 

the 

exterior 

ply. 

average 

buck i i ng 

strength 

ac  t  ua 1  1  y 

increased  slightly.  But,  the 

standard  deviation  of 

buckling 

load  also  increased,  indicating  loss  consistency  in  panel 
strength. 
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3.  Higher  void  contents  also  caused  the  global  panel 
buckling  strength  to  increase.  This  can  be  understood  by 
considering  that.  panel  thickness  increased  linearly  with 
percent  voids,  while  the  buckling  load  is  seen  to  increase 
by  the  square  of  panel  thickness  for  this  geometry.  Higher 
void  contents  are  not  desirable,  however,  as  an  increase  in 
the-  void  content  of  a  panel  increases  the  chance  for 
moisture  absorption  and  subsequent  material  degradation. 

4.  be  1  am  1 na 1 1 ons  in  a  composite  panel,  have  the  potential 
for  more  catastrophic  failure  than  surface  damage. 


De lam l nat i ons 

are 

usual  1 y 

not  visually  observable 

and 

can 

grow  through 

the 

lam  mate 

with  cyclic  loadings. 

For 

the 

de lamination 

sizes 

s  t  ud i ed 

in  this  thesis,  only 

si  i  ght 

decreases  in  global  panel  buckling  were  noted.  However,  in 
general,  impact  on  a  composite  panel  can  potentially  cause  a 
number  of  de 1  am l na t i ons  through  the  pane)  thickness.  This 
would  decrease  panel  strength  even  further  and  would 
increase  the  severity  of  delamination  growth.  Buckled 
sublaminat.es  also  have  relatively  large  deformations  in  the 
z  direction,  which  could  seriously  affect  the  aerodynamics 
of  the  structure,  and  could  conceivably  result  in  aircraft 
instabilities. 

5.  The  phenomenon  of  local  sub  1  am  l  na  t.  e  snapping  of  the 
delaminated  plies  can  be  predicted  using  superposition  ot  a 


r. 


1  1  3 


series 


of"  linear  t"  i  n  1  1 1*  element  models. 


The  use  of  these 


models  results  in  a  fairly  accurate  and  inexpensive 
prediction  of  local  instabilities  caused  by  local  buckling. 
It  is  thought  that  this  analytic  technique  could  be  easily 
extended  to  shells  with  different  geometries  and  loading 
I'undi  *  ions  t  .  >  get  a  bail  park  estimate  of  shell  failure 
1 oads  . 

6.  Seating  loads  should  be  applied  to  panels  prior  to 
experimental  testing.  These  loads  should  be  applied  before 
all  fixture  supports  are  tightened  to  decrease  the  chance  of 
non - s ymme t r 1 c  loading  and  moments  induced  by  differential 
fixture  movements.  Care  should  be  taken  to  ensure  that  the 
panel  is  trimmed  square  and  that  the  load  is  evenly 
introduced  so  that  non - s ymme t r l c  buckling  can  be  avoided. 

7.  The  STAGSC- 1  bifurcation  values  for  the  curved 
cylindrical  panels  could  be  predicted  by  using  a  regression 
curve  that  related  panel  thickness  to  panel  buckling  loads. 
Relations  were  also  found  to  exist  for  other  geometries  such 


as  a  flat  plate 

and 

a  c  i 

rrular 

disc.  1 t 

is  thought 

that 

similar  relations 

cou  1  d 

be 

eas  l 

1  y  est.abl 

l s  h  e  d  for 

other 

shell  geometries, 

and 

that 

the 

amount  o 

f  computer 

t  i  me 

required  to  predict  buckling  of  a  various  thicknesses  for  a 
given  shell  geometry  could  be  greatly  reduced. 


8.  'IV  t' Ion  'lists  proved  to  he  an  e\ee|  lent  Oehurnl  i  ra’ 

material  fur  creating  delaminat  ions,  with  1  U  0  perrent 
delaminat.ion  observed  in  ail  samples.  However,  when  the  0.5 
mil  teflon  was  plated  adjacent  to  the  exterior  ply,  a 
blister  effect  occurred.  This  blister  was  easily  visible  on 
the  surface  of  the  panel ,  and  contained  fractures  in  the 
matrix  of  the  delaminat ion  running  parallel  to  the  liber 
direction.  Initial  finite  element  models  were  run  to 

determine  the  possible  cause  of  this  phenomenon.  There  was 
no  indication  in  the  model  results  that  a  temperature 
differential  during  curing  caused  the  ply  failure.  It  is 
postulated  that  volatiles  formed  during  panel  curing  were 
trapped  beneath  the  teflon  disc  resulting  in  a  vapor 
pressure  which  caused  an  air  bubble  to  form.  A  more 
thorough  investigation  of  this  problem  needs  to  be 
a  1. 1  emp  t  ed  . 

d .  Normal  stresses  at  a  delamination  were  calculated 

through  the  (0/-45)  cross  section  of  a  STAGSC- 1  circular 
plate  model.  These  normal  stresses  were  then  compared  to  an 
average  stress  through  a  (0/-45)  cross  section  of  a  plane 
strain  model  subjected  to  the  same  value  of  strain.  The 
stresses  were  within  0.88  percent  of  each  other.  This 
indicates  that  the  assumption  that  the  delaminated  region  of 
a  curved  panel  can  be  modeled  using  a  plane  strain  model  is 
a  valid  approximation. 
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the  /.  -d  i  ri'ct  i  on  at  t  he  <  •  rark  t  i  p  hth  ea  1  «-u  l  at  e*d  .  Thrsr 
stresses,  normal  to  the*  plane*  of  t.he*  del  ami  nut  ;,  n,  raruied 
from  0 . rt  H  to  1.98  ps l  .  This  level  of  stress  is  much  less 
than  the  stress  required  to  initiate?  crack  drouth  in  this 
material . 


APPENDIX  A 


Constitutive  Matrices  and  Test  Data 


The  Point  Stress  Laminate  Analysis  code  ,  sij:>  111],  was 
used  to  determine  the  stiffness  matrices  for  AS  4  /  3  5  o  1  -  o 
graphite  epoxy.  For  the  quas  1  -  1  so  t.  rop  1  c  symmetric 

t  0 / -  4  5 /  4  5 / 9  0  )  layup,  the  coupling  stiffness  matrix  B  ^  is 
equal  to  zero.  The  matrices  are  computed  for  a  laminate 
thickness  of  0.0384  inches  (corresponding  to  a  ply  thickness 
of  0.0048  inches)  so  that  a  comparison  can  be  made  to  results 
reported  in  previous  work  with  cylindrical  composite  panels 
[7,15].  Matrices  are  then  computed  for  the  median  thickness 
(0.0416)  of  the  panels  tested  in  this  study.  Finally  values 
are  shown  for  a  (0/-45)  layup  so  that  a  plane  strain 
evaluation  could  be  made  of  the  delaminated  region  of  the 
pane  1 . 


For  h=0.0384  inches,  quasi-isotropic  (ply  thickness  t=0.0048) 


315.730  92.535 


0 . 000 


92.535  315.730  0.000  X  10  lb/in 


0 . 000 


0 . 000  111.  600 


63.765 


9.605  -3.867 


9.605  17.362  -3.867 


1  b-  i  n 


-3.86 


-3.86' 


11.947 


F  • 1  r  h  ^  0 . 1  <  4  1  b  ,  (  0  /  -  4  5  /  4  5  /  ‘  K)  )  P  l  v  thickness  t  -  0 . 0  0  o 

s 


X 

l  .) 


3  4  2.043  10  0.247 

100.247  342.043 

0.000  O.UUO 


0 . 0  0  0 
0 . 000 
1  2  0 . 8  y  8 


X  1  0 


B  =  0 

i  0 


D  . 
i  ) 


8  1.071 
12.211 
-4.016 


12.211 
22.074 
-4.9  16 


-4.9  16 
-4.916 
15.190 


1  b-  l  n 


For  (0/-45),  Ply  thickness  t=0.0052  in: 


1 30 . 966 
25 . 062 
-22.728 


25.062 
40 .056 
-22 .728 


-22 .728 
-22.728 
30 . 225 


X  103 


B. 
i  J 


'-172.161 

53.978 


53.978  -59.092 
64.206  -59.092 


-59.092  -59.092  53.973 


lb 


■>  ■ 


Ib/in 


i  b/ 1  n 


D.  . 

1  J 


1 . 18044  0 . 22589  -0 . 20485 

0.22589  0.36103  -0.20485 

-0.20485  -0.20485  0.27242 


1  b-  i  n 


For  (0/-45)  ply  thickness  t= 0.00 4 8  in.: 


A  . 
i  J 


120 . 890 
23.134 
-20.979 


23.134 
36 .974 
-20.979 


-20 .979 
20.979 
27 . 900 


X  103  lb/in 


B. 
i  J 


146.690 

45.993 

-50.350 


45 . 993 
5  4.708 
-50 . 350 


-50 . 350 
-50. 350 
45 . 993 


lb 


D 

i  J 


0.92845  0. 17767  -0.16112 

0.  17767  0.28396  -0. 161  12 

-0. 16112  -0.16112  0.21427 


1  b-  i  n 


The  lab  measurements  and  experimental  test  data  for  this 


thesis  are  recorded  in  the  tables  on  the  following  pages. 


PANEL 

t  (mis) 

%  VOIDS 

EXPER  ( lbs) 

STAGSC-1  X 

DIFF 

Q0-00-1 

40.01 

0.0972 

4795 

5263 

8.9 

qo-oo-2 

40.01 

0.0972 

4445 

5263 

15.5 

Q0-00-3 

41.50 

0.4341 

4779 

5662 

15.6 

Q0-00-4 

41.50 

0.4341 

4995 

5662 

11.8 

Q0-00-5 

41.50 

0.9621 

4945 

5662 

12.7 

Q2-12-1 

42.07 

4.3671 

5085 

5819 

12.6 

Q2-12-2 

41.38 

4.3583 

5214 

5630 

7.4 

Q2-12-3 

41.38 

4.3583 

5321 

5630 

5 . 5 

Q2-12-4 

40.76 

2.2760 

4365 

5462 

20.  1 

Q2-12-5 

40.  15 

2.1762 

5100 

5300 

3.8 

Q2-12-6 

40.  15 

2  1762 

5149 

5300 

2.8 

Q2-12-7 

41.61 

1.2224 

5356 

5692 

5.9 

Q2-12-8 

41.61 

1.2224 

5177 

5692 

9.0 

Q2-23-1 

40.69 

0.1666 

4692 

5443 

13.8 

Q2-23-2 

40.69 

0. 1666 

4570 

5443 

16.0 

Q2-23-3* 

40.80 

-0.0854 

4697 

5473 

14.2 

Q2-23-4* 

40.80 

-0.0854 

4738 

5473 

13.4 

Q2-23-5* 

41.90 

2.0222 

4703 

5772 

18.5 

Q2-23-6** 

41.62 

0.3211 

5325 

5695 

6 . 5 

Q2-23-7** 

42.00 

0.4244 

5160 

5800 

11.0 

Q2-67-1 

40.20 

0.3775 

4233 

5313 

20.3 

Q2-67-2 

40 . 20 

0.3775 

4673 

5313 

12.0 

Q2-67-3* 

41.10 

0.4962 

5357 

5554 

3 . 5 

Q2-67-4* 

41.10 

0.4962 

4255 

5554 

23 . 4 

Q2-67-5** 

42.00 

0.4244 

5295 

5800 

8.7 

Q2-67-6** 

41.62 

0.3211 

5130 

5695 

9.9 

Table  A. 1  -  Panel  Test  Data  (continued  on  next  page ) 


PANEL  t  (nils)  X  VOIDS  EXPBR  (ibe)  STAGSC-1  X  DLFF 


Q2-78-1 

42.61 

3.5801 

d567s 

5969 

6 .  / 

Q2-78-2 

42.61 

3.5801 

5524 

5968 

7 . 5 

Q2-78-3 

41.84 

0.8887 

4795 

5755 

16.7 

Q2-78-4 

41.84 

0.8887 

5659 

5755 

1.7 

Q4-12-1 

41.90 

4.4354 

5229 

5772 

9.4 

Q4-12-2 

41.90 

4.4354 

4885 

5772 

15.4 

Q4-23-1 

39.30 

0.4781 

5224 

5078 

-2.9 

Q4-23-2 

39.30 

0.4781 

4823 

5078 

5.0 

Q4-67-1 

39.60 

0.7278 

4437 

5 1 56 

14.2 

Q4-67-2 

39.60 

0.7278 

4500 

5156 

11.4 

Q4-78-1 

41.77 

0.1169 

4787 

5736 

16.5 

Q4-78-2 

41.77 

0.1169 

4851 

5736 

15.4 

NOTES:  # 

-  Nonsymetric 

buckling  mode 

*  • 

-  3000  lb  seating  load 

t* 

-  2  pieces  of 

teflon,  3000  lb 

seating  load 

For 

SN  1-12  SGI  and  SG2  at  bottom 

of  panel 

Table  A.l  -  Experimental  Data  continued 
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Appendix  B 

Sample  Computer  Input  Decks 


A  researcher  can  invest  an  extremely  large  amount  of 
time  learning  a  computer  operating  system,  and  then  learning 
the  various  computer  codes  necessary  to  model  his  problem. 
This  section  of  the  thesis  is  presented  as  a  guide  for  those 
not  familiar  with  STAGSC-1  or  the  NOS  operating  system  found 
on  the  AFWAL  Cyber  computer  system.  A  sampling  of  STAGSC-1 
input  data  decks  for  various  shell  geometries  and  SQ5  input 
decks  are  presented  and  briefly  discussed  along  with  some 
common  commands  that  a  NOS  user  will  find  useful. 

Most  of  the  jobs  that  were  run  on  the  Cyber  for  this 
research  were  STAGSC-1  finite  element  models.  Since  a 
typical  linear  STAGSC-1  run  for  this  thesis  used  50-150  cp 
seconds  of  computer  time,  these  jobs  must  be  submitted  as 
batch  jobs.  The  following  batch  job  (fig  B.l)  can  be  given 
a  filename  such  as  bjob,  and  then  submitted  using  the 


term i n  a  1  . 


/JOB 


WILDER.P2. 


/USER 


CHARGE , * . 


SETTL ( * ) . 


ATTACH , STAGS 1/UN=D820090. 


ATTACH , STAGS 2 /UN  =  D820090 . 


GET , f i 1 ename. 
STAGS  1 , f i 1 ename . 


RETURN , STAGS  1 . 


STAGS2 . 


RETURN, STAGS2. 

REWIND , * . 

ROUTE , OUTPUT , DC=PR , UN=AL , U JN=BW I LDER , ST=CSA . 


/EOR 


Figure  B.l  -  Typical  STAGSC-1  Batch  Job 


The  batch  job  shown  above  should  be  changed  for  another 
user.  The  new  user’s  nome  should  replace  WILDER  and 
8WILDER,  and  the  filename  of  the  STAGSC- 1  input  deck  to  be 
run  should  replace  filename.  The  batch  job  shown  is  set  up 
to  send  the  output  to  the  printer.  Alternatively,  the 
output  could  be  sent  to  a  file  by  changing  the  ROUTE 
command . 


The  STAGSC- 1  input  deck  shown  below  is  the  input  that 
was  used  to  analyze  the  linear  buckling  load  for  the  curved 
cylindrical  panels  used  in  this  study.  This  model  is  a 
model  for  an  eight  ply  laminate  with  0.0048  inch  thick 
plies.  The  {$(  seen  in  the  input  deck  is  a  divider  between 
fortran  input  and  user  comments.  The  $C1  card  is  the  load 
multiplier  for  the  system  and  multiplies  the  load  entered  on 
the  SQ3  card  to  compute  the  total  system  load. 


LINEAR  ANALYSIS  -  12  IN  PANEL  -  BIFURCATION 
1,1, 1,1, 0,0,1 

1 

1,0,1 

1.0 

1,0,650,0,0 
1,1.0 
19,  19 
1 

18. 844E06, 0.02 18128, .91E06, 1 . , 1 . , 1 . 468E06 , 1 . 

1,1,8 

1,0.00480,0.0 
1,0.00480,-45.0 
1,0.00480,45.0 
1,0.00480,90.0 
1,0.00480,90.0 
1  ,0.00480,45.0 
1 ,0.00480,-45.0 
1 ,0.00480,0.0 
5 

0.0,12.0,0.0,57.3,12.0 

1 

410 

0,0 

100,000 
1 10, 100 
000 , 000 
110, 100 
1 

1,1,0 

0.001,-1,1 ,1,0,0 

2, 2, 2, 0,0, 2 


DELTAS. 001  -  Cl=l 
$B  1  ANALYSIS  TYPE 
SB2  9  SHE  LI.  UN  LI'S 
$B3 

$C1  LOAD  1NCR 

$D2  EIGENVALUE  CNTRL 

$D3  CLUSTER  DEFN 

$F1  NO.  ROWS  COLS 

$11  MATL  #1 

$12  MATL  PROP 

$K1  8  LAMINAE 

$K2  PLY,  T,  ORIENT 

SK2 

$K2 

$K2 

$K2 

$K2 

$K2 

$K2 

$M1  SHELL  TYPE 
$M2A  SHELL  GEOMETRY 
$M5 

$N  1  ELEMENT  TYPE 
$P1 

$P2  BOliNDARY  OOND 

$P2 

$P2 

$P2 

$Q1 

$Q2 

$«3  LOAD 

$R1  aTR’l  FORMA! 


Figure  B.2  -  STAGSC- 1  Input  Dock  for  a  Curved  12  inch  Panel 


-  v.---.  .■.'-'.■.••.•vyvyyu.'uu-fv-ucu;.-.; 


The  input  deck  shown  in  figure  B.3  below  was  one  of  the 
input  decks  used  in  this  study  to  perform  t.he  convergence 
study  for  buckling  ot'  a  circular  plate  under  a  uniform  axial 
displacement.  This  particular  input  deck  is  for  a  10x8,  90 
degree  model.  The  material  property  input  card  $11!  shows 
the  input  required  for  an  isotropic  material.  The  shell 
geometry  card  $M1  is  the  card  which  defines  this  shell  as  an 
annular  ring.  The  $ M 2 A  card  gives  the  radius  of  the  inner 
boundary,  the  radius  of  the  outer  boundary,  and  the  angles 
in  degrees  swept  out  by  the  two  side  boundaries.  The 
element  type  chosen  for  this  model  was  the  QUAF  410 
element.  Using  this  element  and  the  given  plate  geometry, 
STAGSC- 1  will  automatically  generate  each  element  in  the 
mesh  at  intervals  defined  by  the  user.  The  four  $P2  cards 
are  the  cards  which  define  the  constraints  at  the  shell’s 
boundaries  in  the  following  order:  displacement  in  the  x,v,z 
directions,  rotations  with  respect  to  x,  y,  z.  A  one 
signifies  freedom  to  move  while  a  zero  signifies  a 
constrained  node.  The  user  has  the  option  of  specifying  a 
separate  load  and  constraint  at  each  node  as  can  be  seen  in 
figure  B.5  or  of  maintaining  the  same  constraints  and  loads 
at  every  node  on  a  boundary  as  is  done  in  figure  B.4.  For 
this  ease,  the  third  card  specifies  a  clamped  edge  condition 
for  the  outside  arc  of  the  model.  The  displacement  load  is 
deft  iied  by  the  $W3  card.  This  card  not,  only  gives  the 
d i s p 1 u< emen t  value  but  the  boundary  to  which  it  applied. 


c 


LINEAR  ANALYSIS  -  2IN  DISK  OPLY 
1,1, 1, 1,0,0,  1 
1 

1.0,1 

1.0 

1,0,650,0,0 
1, .001 
10,8 
1 

18. 844E06 ,0.3 
1,1,1 

1,0.00480,0.0 

4 

00.0,1.0,0.0,90.0 

1 

410 

0,0 

001,111 

101,011 

110,101 

101,011 

1 

1,1,0 

-.001,-1,1,10,0,0 
1,1, 1,0, 0,1 


-  BIFUR  -  DISK5  -  ISOTROPIC 

$B1  ANALYSIS  TYRE 
$B2  NO.  SHELL  UNITS 
$B3 

SCI  LOAD  I NCR 

SDL  EIGENVALUE  CNTRL 

$D3  CLUSTER  DEFN 

$F1  NO.  ROWS  COLS 

$11  MATL  #1 

$12  MATL  PROP 

SKI  8  LAMINAE 

$K2  PLY  THICK,  ORIENT 

SMI  SHELL  TYPE 

$M2A  SHELL  GEOMETRY 

$M5 

$N1  ELEMENT  TYRE 
$P1 

$P2  BOUNDARY  OOND 

$P2 

$P2 

$P2 

$Q1 

$Q2 

$Q3 

SRI  OUTPUT  FORMAT 


Figure  B.3  -  STAGSC-1  Input  for  Two  Inch  Circular  Plate 


The  next  input  deck  seen  in  figure  B.4  below  was  used 
to  generate  a  STAGSC-1  circular  plate  model  similar  to  the 
model  used  by  Shivakumar  and  Whitcomb  [13]  in  their  work 
with  elliptical  de 1 am i na t i ons  (see  fig  5.7).  This  model  was 
more  difficult  to  generate,  and  didn’t  agree  as  well  with 
the  closed  form  solution  as  the  10x8  model.  A  couple  of 
features  of  this  model  to  note  are  the  input  cards  $N9A  and 
$N9B.  These  cards  generate  a  series  of  similarly  oriented 
elements  around  the  model.  The  first  card  defines  the 
number  of  cards  necessary  to  define  the  model,  and  the  other 
cards  define  the  element,  the  number  of  elements  in  the 


series  and  t  h<-  node  numbers  which  define  the  first  element 
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LINEAR  .ANALYSIS  -  2  IN  DISK  OPLY 
1,1 ,1,1 ,0,0,1 
1 

1,0,1 

1.0 

1,0,650,0,0 

1,0.001 

9,33 

1 

18. 844E06, 0.02 18128, .91E06, 1. , 1. 

1,1,2 

1,0.00480,0.0 

1,0.00480,-45.0 

4 

0.0,1.0,0.0,360.0 

1 

0,0,0 

15 

320.1.8.1.1.2.1.2.5 

320.1.8.2.1.3.3.2.5 
320,1,8,2,1,3,1,3,3,0,4 
320,1,8,2,5,3,3,3,5,0,4 
320,1,16,3,1,4,1,4,2,0,2 
320,1,16,3,1,4,2,3,3,0,2 
320,1,16,4,2,4,3,3,3,0,2 

320.3.16.4.1.5.1.4.2.2.2 

320.3.16.5.1.5.2.4.2.2.2 

320.3. 16.4.2.5.2.5.3.2.2 

320.3.16.4.2.5.3.4.3.2.2 

320.2.16.5.1.6.1.6.2.2.2 

320.2.16.5.1.6.2.5.2.2.2 

320.2.16.5.2.6.2.5.3.2.2 

320.2.16.6.2.6.3.5.3.2.2 

0,0 

001,111 

101,011 

110,010 

101,011 

1 

1,1,0 

-.001,-1,1,9,0,0 

2, 2, 2, 0,0, 2 


BIFUR  -  DISKA  -  0/-45  PLIES 

$B1  No.  SHE  Li.  UNITS 
$B2  NO.  SHELL  IN  ITS 
$B3 

$C1  LOAD  INCH 
$D2  EIGENVALUE  U'NTKL 
$D3  CLUSTER  DEFN 
$F1  NO.  ROWS  COLS 
$11  MATL  #1 

1 . 468E06 , 1 .  $12  MATL  PROP 
$K1  #  LAMINAE 
$K2  PLY  THICK,  ORIENT 
$K2 

$M1  SHELL  TYPE 
$M2A  SHELL  GEOMETRY 
$M5 

$N1  ELEMENT  TY  PE 

$N9A 

$N'9B 

$N9B 

$N9B 

$N9B  ELEMENT!'  ORIENTATION 


$P1 

$P2  BOUNDARY  COND 

$P2 

$P2 

$P2 

$Q1 

$Q2 

$Q3 

SRI  OUTPUT  FORMAT 


« 


Figure  B.4  -  STAGSC-1  360  Degree  Circular  Plate  Mode 1  Using 
the  TRINC  320  Element. 
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1  28 


method  of  loading  individual  nodes  on  a  situ,  fur--  urn 
■'.'onstra  in  1  ng  others.  Note  that  two  load  h  yst  eras  hi--  H  in 
by  the  vise  of  two  $\jz  cards  f  .>  1  lowed  h\  a  series  .  >  f  *w-:  !  ■ 
records  . 


LINEAR  ANALYSIS  -  21N  DISK  0P1.Y  -  B 1  FUR  -  D 1  SKA  - 

U/-45  PLIES 

1,1, 1,1, 0,0,1 

$B1 

NO.  SHK.1J,  l  MIN 

1 

$B2 

NO .  SHE11.  1  N I  I'M 

1,0,1 

SB3 

1 .0 

SCI 

LOAD  INCH 

1 ,0,650,0,0 

SD2 

EIGENVALUE  ( INTEL 

1 ,0.001 

$D3 

CLUSTER  DEKN 

y  ,33 

SKI 

NO.  ROWS  COLS 

1 

$11 

MATL  ttl 

18 . 844E06 ,0.0218128, . 8 1E06 , 1 . , 1 . , 1 .468E06, 1 . 

$12 

MATL  PROP 

1,1,3 

SKI 

#  LAMINAE 

1,0.00480,0.0 

SK2 

PLY  THICK,  OklEN’l 

1,0.00480,-45.0 

SK2 

4 

SMI 

SHELL  TA  PE 

0.0,1.0,0.0,360.0 

SM2A  SHELL  UEUMETRY 

1 

SM5 

0,0,0 

$N1 

ELEMENT  TYPE 

15 

SN9A 

320, 1,8, 1,1, 2, 1,2, 5 

SN9B 

320,1,8,2,1,3,3,2,5 

SN9B 

320,1,8,2,1,3,1,3,3,0,4 

SN9B 

320,1,8,2,5,3,3,3,5,0,4 

SN9B  ELEMENT  ORIENTATION 

320, 1,16,3,1,4, 1,4, 2, 0,2 
320,1,16,3,1,4,2,3,3,0,2 
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Figure  B.5  -  STAGSC- 1  Model  of  Circular  Plate  With  Multiple 


Load  Systems  (Continued  on  Next  Page ) 


FikTure  B .  5  Continued 

$10 

'('1,11! 

$P2  But  NO.AKV  O  >Nl) 

101  ,'M  1 

$P2 

110.101 

$P2 

101,011 

$P2 

- 

$Wl 

i ,  .<3 ,  j 

$0)2 

-,0Ol ,-l 

1 ,0,8,0 

SQ.i 

-.001 , - ! 

1 ,0,0,0 

$W3 

-.001  ,-l 

1,0,  i  U  ,  0 

$0)3 

0 . 000 , -  1 

1  ,  9  ,  1 . 0 

$0)3 

0 . 000, -  1 

1 ,0,2,0 

SW3 

0 . 000 , -  1 

1 ,0,3,0 

$0)3 

0 .000 , -  1 

1 . 0 , 4 , 0 

$0)3 

0 . 000 , -  1 

1 , 0 , 5 , 0 

$0)3 

0.000,-1 

1 , 0 , 6 , 0 

$0)3 

0 . OOO , -  1 

1  ,9, 7 ,0 

$0)3 

0.000,-1 ,1,0,11,0 

$W3 

0 . 000 , -  1 

1,9,12,0 

$y3 

0.000,-1 

1,9,13,0 

$W3 

o.ooo.-i 

1.9,  14,0 

$0)3 

0 . 000 , -  1 , 

1,9,15,0 

$W3 

0.000,-1 

1,9,16,0 

$0)3 

0 . 000 , -  1 , 

1,9,17,0 

$W3 

0.000,-1 , 

1  ,  9  ,  1  H  ,  0 

$Q3 

o.ooo.-i , 

1,9,19,0 

$0)3 

0.000,-1 , 

1 ,9.20,0 

$0)3 

0 . 000 , -  1  , 

1,9,21,0 

$0)3 

0 . 000 , -  1 , 

1,9,22,0 

$o)3 

0 . 000 , -  1  , 

1.9,23,0 

sy3 

0 . 000 , -  1  , 

1 ,9,24,0 

$Q3 

0 . 000 , -  1  , 

1,9,25,0 

$0)3 

0.000,-1, 

1,9,26,0 

$0)3 

0.000,-1  , 

1,9,27,0 

$0)3 

0 . 000 , -  1  , 

1,9,28,0 

$0)3 

0.000,-1 , 

1,9,29,0 

$0)3 

0 . 000 , -  1  , 

1,9,30, 0 

$0)3 

0.000,-1  , 

1,9,31,0 

$0>3 

0.000,-1  , 

1,9,32,0 

$0)3 

0.000,-1 , 

1,9, 33 , 0 

$0)3 

2,  1  ,0 

$0)2 

o . ,  - 1 , 2 ,  y 

,0,0 

$0)3 

2,2,11,0,0 

') 

t 

SR  1  <XTMT  H>KMA1 

FiKurt' 

B.5a  (Continued)  -  STAGS*- 1 

lri'uUr  I’latc  *  ;  •  h 

Mu i i i p 1 e 

Load  S  vst  » > m s 

1  .io 

v  ^.v. -y-  .-03.  [v '.  oN ' .  ; A-i 

The  Point,  Stress  Laminate  program  Sy5  was  useu  t.o 
determine  the  constitutive  relations  used  in  this  research. 
A  sample  input  deck  is  given  in  the  figure  below.  Sw5  can 
be  used  to  compute  most  of  the  values  necessary  t.o  perform 
classical  laminate  calculations.  Models  can  be  run  which 
compute  force  resultants  and  moment  resultants  for  a 
temperature  differential,  which  compute  resultants  for  a 
given  strain,  and  which  compute  average  material  properties 
for  a  given  laminate. 

PROEERT 1 KS  Oh  A  LAMINATE  <0,-+45,90)S 
C  0  0  0  0  2  10 

18844000.  1468000.  0.28  0 10000.  0.  0. 

1  1  0.  0.0048 

1  -45.  0.0048 

Figure  B.6  -  Sample  SQ5  Input  Deck 

,  select  ion  ot  common  NOS  commands  used  on  the  AFWAL 
■  r>>  r  m futer  system  is  presented  here  for  subsequent 
-O  'liter:  *  w  n  >  must  learn  to  use  STAGSC-i  on  the  current 
s  -n  .  M  is  first.  suggested  that  the  new  user  load  the 

’■  1  I  1  1  r  •  f  i  I  •  inti  h  i  s  accoun  t  by  fn  l  i .  wing  the  d  i  re«' t  ions 


outlined  in  the  "An  introduct  1  on  to  AK1T  ('omputer  Systems 
pamphlet..  The  Proof  lie  contains  a  library  of  abbreviated 
commands  which  are  similar  to  UNIX  c< immands .  A  help  file  is 
also  available  which  gives  a  summary  ot  most  of  the 
abbreviated  commands  and  which  can  be  sent  to  the  printer  so 
that  the  user  has  a  handy  reference. 

get ,  filename  -  retrieves  a  file  from  the  library  into 
the  working  environment.  The  file  must  be  retrieved  into 
the  working  environment  before  it  can  be  operated  on  or 
executed  by  the  user. 

save  ,  f  i  1  ename-  new  f  i  1  ename  -  used  after  an  editing 
session  to  save  the  new  version  of  a  file  to  the  user's 
1 ibrary 

replace , f i 1 ename  -  used  after  an  editing  session  to 
replace  an  old  file  with  a  new  version  in  the  library 

rts , filename  -  routes  a  file  to  the  printer 

submit ,  batch  job, to  -  submits  a  batchjob  and  sends  the 
dayfile  to  the  waitque. 


1  A  2 


0 


gives  the  status  ot'  all  current  executing  jobs  and 


t'  L 

tells  which  .jobs  are  in  the  wait  queue. 

qg  ,  .jobnarae  -  retrieves  a  job’s  day  file  from  the 
uaitque.  The  jobname  is  a  four  character  code  assigned  by 
the  computer  which  defines  the  job  and  is  seen  using  f 1 . 

pur  ate,  fi  lennme  -  deletes  a  file  from  the  directory 

Is  -  a  Proof  ile  command  which  shows  all  ot'  the  files  in 
the  user’s  directory 

mor e,  filename  -  prints  a  file  to  the  screen  one  page  at 
a  time.  Hitting  the  space  bar  or  the  carriage  return  bring 
up  the  next  page  of  the  file. 
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physically  modeled  in  the  panels  by  removing  a  portion  of  the 
exterior  plies.  The  panels  were  then  tested  by  compress i ve 1 y 
loading  them  and  a  comparison  was  made  to  buckling  predictions 
obtained  using  a  STAGSC-1  shell  program.  These  tests  indicated 
that  panels  which  have  suffered  minor  surface  damage  do  not 
deviate  significantly  from  buck  1 i ng  predictions  obtained  using  a 
STAGSC-1  linear  bifurcation  model. 

Panels  were  also  tested  which  had  varying  thicknesses  and 
variations  in  void  content  due  to  faulty  manufacturing 
techniques.  These  panels  were  also  compressively  loaded,  and  it 
was  found  that  high  void  content  increased  panel  thickness,  which 
resulted  in  higher  compressive  strengths. 

Composite  laminates  subjected  to  a  low  speed  impact,  such  as 
a  dropped  tool  or  a  manufacturing  load,  often  develop  an  internal 
delamination.  This  delamination  may  result  in  the  reduction  of 
the  panel's  strength  when  subjected  to  compressive  load. 

Since  curved  panels  are  3-dimensional,  and  buckling  is  a  non¬ 
linear  phenomenon,  the  compressive  load  which  will  cause  curved 

panels  to  become  unstable  is  extremely  hard  to  predict 
analytically.  This  thesis  presents  a  technique  whereby  the  local 
buckling  loads  at  the  delamination  may  be  predicted  using  a  2- 

dimenaional  model  with  a  plane  strain  correction.  This  model 

yielded  predictions  of  local  instability  within  30  percent  of 

experimental  values. 


UNCLASSIFIED 


